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vortex-lattice  simulation.  The  spanwise  jet  is  shown  to  augment  the  stability 
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Certain  scale  or  Reynolds  Number  effects  are  formulated,  particularly  concerning 
the  prediction  of  flow  stability  such  as  vortex  burst  over  the  wing  or  in  the 
near  wake.  The  preliminary  results  indicate  that,  possibly,  such  instabilities 
are  less  likely  to  occur  for  vortex  lift  control  in  full-scale  flight.  In  other 
respects  the  differences  between  model  ana  large  scale  theoretical  predictions 
are  found  to  be  negligible,  as  is  generally  confirmed  by  the  high  Reynolds  Number 
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FOREWORD 

This  research  on  a  unique  lift-augmentation  concept  was  performed 
under  Contract  No.  N00014-72-C-0237  for  the  Office  of  Naval 
Research.  Mr.  T.  L.  Wilson  was  Technical  Contract  Monitor  for 
Aeronautics,  Code  461.  The  results  are  reported  in  the  following 
volumes. 

VOLUME  I  -  EXPERIMENTAL  RESEARCH 

VOLUME  II  -  THREE-DIMENSIONAL  THEORY  FOR  VORTEX-LIFT 
AUGMENTATION 


SUMMARY 


An  attempt  is  made  in  this  volume  to  develop  and  demonstrate  an  analyti¬ 
cal  procedure  by  which  the  phenomena  of  vortex-1 1ft  control,  using  spanwise 
blowing  or  mass/momentum  addition,  can  be  predicted  to  a  useful  degree  of 
accuracy.  The  approach  is  predicated  on  the  assumption  that  the  simplest 
conceptual  and  heuristic  reasoning  is  to  be  used  where  alternatives  of  more 
complexity  may  possess  equal  or  greater  plausibility.  Within  this  context 
a  measure  of  success  was  achieved  in  terms  of  understanding  the  physical, 
viscid  and  inviscid,  mechanisms  of  coaxial  jet-and-vortex  augmentation  and 
entrainment  as  the  pair  are  swept  outboard  and  aft  over  a  lifting  surface. 

As  anticipated,  limitations  ?n  the  theory  became  apparent  in  early  stages; 
but  even  these  results  are  useful  in  pointing  out  the  requirements  for 
future  research,  as  well  as  suggesting  even  broader  vortex  theory  applications 
such  as  possible  techniques  for  controlling  or  modifying  the  growth  and  dis¬ 
sipation  of  potentially  hazardous  wing  tip  vortices. 

The  theoretical  analysis  applies  two-dimensional,  potential  flow 
solutions  using  conformal  mapping  to  obtain  vortex-sink  strengths  as  starting 
conditions  for  a  three-dimensional  vortex  lattice  simulation.  This  latter 
computer  program  is  modified  to  account  for  source/sink  and  viscous  vortex 
interactions  with  a  crossflowing  jet.  Nearly  optimum  conditions  for  vortex- 
lift  augmentation  are  selected  for  emphasis  in  the  simulation.  Thus,  an 
under-expanded  (supersonic)  jet  is  assumed  to  issue,  over  a  wing  at  stall 
attitudes,  from  a  nozzle  pointed  in  an  approximate  spanwise  direction  at  a 
chordwise  position  known  to  allow  a  sizeable  Helmholtz  flow  to  normally 
persist  due  to  leading  edge  separation.  The  boundary  conditions  result  In  a 
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vortex  forming  above  and  aft  of  the  leading  edge,  as  discussed  In  the  experl-  * 

mental  work  (Volume  I);  and  the  three-dimensional  solution  provides  vortlclty 

5 

augmentation  through  a  Stokes1  Integration  of  the  jet-boundary  layer  (wall- 
je'c)  shear  profiles.  Empirical  hypotheses,  such  as  the  Polhamus1  leading 
edge  suction  analogy,  are  Introduced  to  justify  a  rapidly  converging 
iterative  procedure  which  yields  numerical  solutions. 


Definitive  results  are  obtained  on  wing  lift  and  corresponding  pressure 
distributions  for  comparison  with  the  test  data  presented  In  Volume  1. 

Certain  scale  or  Reynolds  Number  effects  are  formulated,  particularly  con¬ 
cerning  the  prediction  of  flow  stability  such  as  vortex  burst  over  the  wing 
or  in  the  near  wake.  The  preliminary  results  Indicate  that,  possibly,  such 
Instabilities  are  less  likely  to  occur  for  vortex  lift  control  in  full-scale 
flight.  In  other  respects  the  differences  between  model  and  large  scale 
theoretical  predictions  are  found  to  be  negligible,  as  is  generally  confirmed 
by  the  high  Reynolds  Number  tests  (Volume  l). 
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The  search  for  a  greater  understanding  of  vortex  flows  very  likely 
provides  one  of  the  most  persistent  historical  goals  of  fluid  mechanics 
research.  The  elusive  character  of  vortices,  forming  and  dispersing  in 
sometimes  perturbing  and  often  beneficial  conditions,  is  largely  dependent 
upon  the  stabilizing  or  destabilizing  characteristics  of  the  local  flow 
environment.  Basic  analytical  and  experimental  research  on  the  mechanisms  of 
vortex  flows  in  bluff  body  wakes,  as  described  in  Reference  1,  have  resulted 
in  specific  means  for  controlling  vortex  growth  and  stability.  It  is 
demonstrated  therein  that  the  Strouhal  frequency  in  the  Von  Karman  street  can 
be  predicted  theoretically  if,  and  only  if,  the  three-dimensional  viscous 
Navier-Stokes  equations  are  utilized  so  that  axial  vortex  flow  can  be  con¬ 
sidered.  Various  means  were  used  in  tests  (Ref.  l)  for  modifying  the  axial 
mass-flow  and  vorticity  propagation  characteristics  for  vortices  forming  in 
near-wake  or  separated  flow  regions.  Such  modifications  resulted  in  elimi¬ 
nation  of  the  periodic  vortex  shedding  and  restoration  of  the  Foppl  vortex 
pair  behind  bluff  bodies  at  supercritical  Reynolds  numbers.  Further  experi¬ 
ments  at  that  cime,  reported  in  References  1  through  5  demonstrated  that 
either  suction  or  blowing  parallel  to  the  axis  of  vortices  in  cavities,  or 
separated  flow  regions,  could  stabilize  and  "lock"  an  augmented  vortex  in 
position  so  as  to  produce  direct  lift  increments.  Subsequently,  many 
practical  applications  have  been  found  and  published  (Refs.  2,  6-8)  concern¬ 
ing  both  lift  and  drag  as  well  as  stability  and  control  implications. 

The  present  research  program  extends  the  earlier  research  by  showing 
more  generalized  dependence  of  the  flow  control  mechanisms  upon  the  Rossby 
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number  (ratio  of  1 inear-to-angular  momentum)  and  Radial  Inflow  Reynolds 
Number  previously  discussed  In  Reference  1.  For  example,  the  vortex  growth, 
persistence  and  stability  are  of  obvious  Importance  to  the  spanwlse  blowing, 
vortex  momentum  exchange  phenomenon.  Thus,  stability  criteria  and  effective¬ 
ness  parameters  are  developed  which  allow  preliminary  estimates  to  be  calcu¬ 
lated  for  vortex- lift  augmentation  applied  to  more  general  aircraft 
configurations.  However,  future  research  is  required  for  complete  resolution 
of  this  problem  for  final  design  purposes. 
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2.  TWO-DIMENSIONAL  POTENTIAL  FLOW  ANALYSIS 


2.1  BASIS  FOR  THE  MATHEMATICAL  MODEL 


2.1.1  Experimental  Observations 


The  extensive  smoke  flow  visualization  experiments  on  spanwise  blowing 


have  been  studied  for  a  wide  range  of  flow  conditions.  A  typical  flow 


pattern  is  shown  in  Figure  2  of  Volume  i.  This  photograph  shows  the  flow 


structure  in  a  chordwise  vertical  plane  near  the  wing  root  with  a  jet  blowing 


into  the  region  which  would  be  fully  separated  at  this  incidence  if  no 


blowing  were  present.  A  distinctly  discrete  vortex  is  seen  downstream  of  the 


leading  edge  separation  point  but  forward  of  the  lateral  jet  flow.  Fluid  is 


observed  to  spiral  in  toward  the  vortex  center,  suggesting  the  presence  of  a 


sink  as  viewed  in  two  dimensions.  Whatever  the  complex  three-dimensional 


interaction  may  be  between  the  jet  and  the  stalled  flow  field,  it  seems 


apparent  that  the  jet  causes  a  vortex  to  appear  and  to  maintain  its  existence 


by  inducing  an  axial  flow  in  the  vortex.  The  apparent  sink  strength  is 


sufficient  to  duct  fluid  mass  out  the  wing  span  at  the  same  rate  at  which  it 


is  entrained  from  the  leading  edge  and  freestream.  The  flow  downstream  of  the 


jet  centerline  is  steady  and  appears  to  reattach  somewhere  along  the  chord  aft  of 


the  centerline.  For  some  conditions  there  are  multiple  stagnation  points 


along  the  upper  surface  of  the  chord.  In  every  case  the  flow  near  the 


trailing  edge  is  well  behaved  so  that  the  Kutta-condi t ion  may  be  applied  as  a 


necessary  condition  in  the  potential  flow  model.  However,  these  observations 


of  flow  pattern  apply  primarily  for  angles  of  attack  above  the  stall  angle  of 


the  airfoil  section  being  considered.  Based  on  the  preceding  visual  evidence, 


a  mathematical  model  is  formulated  as  described  below. 
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2.1.2  Mathematical  Formulation 

Since  the  leading  edge  detachment  is  present  in  practical  cases,  a  flat 
plate  approximation  of  the  airfoil  is  considered  to  be  reasonable.  This  leads 
to  a  simpler  mathematical  representation  in  the  complex  plane  while  retaining 
the  essential  vortex  lift  augmentation  effects  demonstrated  for  the  wing  u.ider 
test.  The  other  elements  of  the  flow  are  a  vortex,  a  source/sink  at  the  center 
of  the  vortex,  and  a  source  to  represent  the  jet.  The  vortex  feeding  sheet 
from  the  leading  edge  is  not  represented  in  the  present  model;  this  leads  to  a 
more  tractable  mathematical  solution  and  is  considered  to  be  justified  in  a 
first  attempt  to  assess  the  effects  of  the  more  notable  flow  mechanisms.  Also, 
the  representation  of  the  jet  in  two  dimensions  as  ..  source  requires  careful 
physical  interpretation  and  this  subject  will  be  deferred  for  later  discussion. 

Conformal  mapping  techniques  are  convenient  for  representing  the  two- 
dimensional  flow  field.  The  airfoil  is  taken  siong  the  real  axis  in  the  Z- 
piane  between  (-2,  +2)  as  shown  in  Figure  1.  The  jet  position,  Zj ,  and  its 
source  strength,  mj ,  are  assumed  to  be  known.  To  represent  typical  chordwise 
sections,  this  implies  a  priori  knowledge  of  crossflowing  -  jet  phenomena  as 
discussed  in  the  section  on  three-dimensional  methods.  The  vortex  center  is 
located  at  Zv,  and  it  has  a  circulation  strength,  Kv;  a  source/sink,  mv,  is 
also  located  at  Zv  (whether  a  source  or  a  sink  is  present  will  be  determined 
by  the  sign  in  the  solution).  The  freestream  is  set  at  angle  a,  as  shown  in 
the  figure,  and  it  has  a  v*>ijcity  magnitude,  Vw.  Finally,  there  is  airfoil 
bound  circulation  of  strength  Ka.*  The  Z-plane  is  mapped  into  the  circle 
plane  or  c-plane  by, 

*For  convenience  a  factor  2ir  is  included  in  the  m's  and  K's,  thus  the  source/ 
sink  strength  in  vol  ./sec. /ft.  is  2?im  and  the  circulation  in  ft.2/sec.  is  2ttK. 


Z-  PLANE 


£  -  PLANE 


Figure  1.  Conformal  Mapping  Coordinates  in  Transform  Planes 

for  the  Two-Dimensional  Analysis  of  Vortex-Lift  Augmentation. 
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and  the  complete  potent!?, 1  for  the  ?-plane  Is  found  by  add<ng  the  complex 
potentials  of  each  of  the  elements  (Reference  9)  as  follows: 

w(c)  -  V  (ce5a  +  I  e"la)  +  i  Kv[ln(;-?  )  -  1 n (c  “=  ) 

C  *v 

+  Inc]  -  mj  [In  (C'Cj)  +  Me-  i  )  -  Inc]  -  mv[ln(c-Cv) 

+  ln(c-  1  )  -  Inc]  +  i  Ka  In  C.  (l) 

Cv 

The  bar  denotes  complex  conjugate  and  cv  ar|d  Cj  are  the  vortex  and  jet 
positions  in  the  C~plane  (Figure  1). 

Two  kinematic  conditions  are  available  for  solution  of  the  problem: 

(l)  as  noted  earlier,  the  flow  appears  to  be  attached  in  the  neighborhood  of 
the  trailing  edge,  so  that  the  Kutta  condition  of  potential  flow  theory  may 
be  applied;  (2)  for  the  "actual"  flow  the  vortex  is  stationary  with  respect 
to  the  airfoil;  therefore,  the  potential  flow  muse  have  a  stagnation  point 
at  the  vortex  center  when  the  singularities  at  that  point  are  removed.  The 
first  condition  generates  one  equation  since  the  trailing  edge  of  the  air¬ 
foil  transforms  into  a  point  on  the  real  axis  in  the  c-plane.  In  general 
the  vortex  location  Cv>  will  be  complex  and  will  generate  two  equations. 

Thus,  the  following  results, 
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where,  in  the  last  equation,  the  singular  parts  at  ?  =  cv  have  been 
removed.  The  real  and  imaginary  parts  of  this  equation  are  set  to  equal  to 
zero  so  that,  along  with  equation  (2),  there  are  three  equations  which  may 
be  used  to  determine  Ka,  Ky,  andmy,  with  mj  assumed  to  be  known. 


2.1.3  Limitations  of  the  Model 

The  difficulties  associated  with  this  model  are  that  leading  edge 
separation  Is  not  represented,  with  the  attendant  large  reduction  in  flow 
speed  over  the  upper  surface  of  the  plate;  and  that  a  uniquely  determinable 
position  does  not  exist  for  the  vortex.  Regarding  the  second  of  these  points, 
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if  the  condition,  (refer  to  Figure  1  for  definition  of  variables): 


is  expanded  mathematically  with  sv  =  rve'®v,  the  resulting  equations  involve 
rv  to  the  sixth  power.  This  implies  that  at  least  one  of  the  quantities,  Ka, 
Ky,  or  mv  must  be  quite  sensitive  to  vortex  position,  with  all  other  condi¬ 
tions  unchanged.  This  will  be  discussed  further  in  the  section  on  numerical 
solutions.  It  seems  clear  from  flow  visualization  experiments  that  super¬ 
position  of  the  several  potential  flows  will  not  fully  represent  the  true 
flow  field  due  to  the  detached  flow  condition  at  the  leading  edge.  The 
inclusion  of  a  vortex  feeding  sheet  from  the  leading  edge  to  the  vortex 
would  allow  a  better  approximation  and,  in  addition,  would  provide  a  dynamic 
condition  for  uniquely  determining  the  vortex  position.  In  lieu  of  this 
complicating  feature,  a  modified  condition  for  stagnation  at  the  vortex 
center  is  used.  The  region  where  the  trapped  vortex  is  finally  established 
is  otherwise  a  fully  separated  wake.  Fortunately,  results  are  available  from 
surveys  of  the  velocity  distribution  in  the  near-wake  of  a  stalled  flat 
plate,  for  a  wide  range  of  angle  of  attack  above  stall,  indicating  that  the 
mean  "cavity"  velocity  is  about  .1  (Refs.  6,  10,  and  11).  Thus,  in  place 
freestream  velocity  in  equation  (3)  for  vortex  stagnation,  the  value 

^cav  =  •  ^  (5) 


is  used  herein.  A  similar  reduction  of  upper  surface  velocity  magnitude  is 
computed  more  analytically  in  the  three-dimensional  solution  because  more 
complete  knowledge  of  the  "cavity"  flow  reversal  properties  is  then  available. 


2 . 2  METHODjDP  SOLUTION 


Once  a  mass  flow  rate  is  established  for  the  jet  at  a  given  angle  of 
attack  ana  free  stream  speed,  the  following  procedure  is  used  to  obtain 
solutions.  Since  the  "actual"  location  of  the  vortex  is  not  known,  a 
matrix  of  positions  of  the  vortex  must  be  computed  and  the  resulting 
solutions  for  sink  strength,  vortex  circulation,  and  airfoil  circulation 
are  then  examined  on  physical  grounds  to  determine  the  range  of  most 
reasonable  values.  While  sucn  analyses  are  valuable  in  providing  an 
unaer standing  01  the  nature  of  the  flow  interactions,  they  will  not  yield 
unique  solutions  since  only  one  location  of  the  vortex,  unknown  because 
the  problem  as  posed  is  under-determined,  corresponds  to  reality. 

The  two-dimensional  representation  of  the  jet  as  a  source  is  intuitively 
appealing.  However,  the  methou  for  deriving  an  equivalent  strength  from 
three-dimensional  data  (Volume  I  of  this  report)  is  not  entirely  clear. 
For  tne  results  presented  here,  the  entrainment  rate  from  laser-doppler- 
velocimeter  (LDV)  measurements  was  used  to  compute  tne  rate  of  expansion 
ol  the  jet  along  its  axis.  This  expansion  then  was  equated  to  a  source 
out-flow  in  cubic-feet-per-second,  per  foot  cl  length.  For  0^  =  .46,  the 
experimental  value  chosen  for  the  most  detailed  theoretical  calculations, 
the  jet  source  strength  is  found  to  be  m^  =  2.16.  The  actual  value  used 
for  the  following  computations  was  rounded  off  to  2.0. 


j  zl ±1 e  ■■ :  A 
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2.3  SiU^I^COMHJTE^  RESULTS 

Figure  2  presents  results  for  the  jet  positioned  at  ,25c  chordwise  and  at 
.075c  above  the  airfo.l.  The  vortex  is  placed  at  .10c,  ,15c  and  .20c  chord- 
wise.  The  position  above  the  airfoil  in  percent-of-chord  is  the  abscissa. 

For  this  case  the  outer  free-stream  speed  is  100  fps.  and  a  =  25  degrees. 
Clearly,  K  and  K  are  not  as  sensitive  to  chordwise  vortex  position  as  is 

cv  V 

the  vortex  sink  strength  m^.  No  criteria  exist  in  the  model  to  Belect  the 
correct  vortex  position  in  this  solution  array;  however,  an  argument  can  be 
made  for  bounding  values  on  physical  grounds.  The  vertex  cannot  have  source 
flow  at  the  core  since  this  is  an  unstable  (or  non-sustainable)  condition 
(Ref.  l).  Therefore,  a  lov/er  bound  on  vertical  vortex  position  is  indicated 
in  Figure  2  for  each  chordwise  position.  The  lower  dotted  line  denotes  the 
circulation  strengtn  at  the  given  angle  of  attack  for  a  fully  stalled  flat 
plate  (Refs.  6  and  12),  and  it  seems  reasonable  to  assume  that  K  will  not  oe 
less  than  this  value.  The  upper  dotted  line  denotes  the  potential  flow 
circulation  limit  for  the  airfoil  alone.  The  lines  above  this  level  are  the 
sums,  +  Ky,  for  the  three  chordwise  position  solutions.  It  is  interesting 
to  note  that  circulation  augmentation  declines  to  zero  for  about  the  same 
vertical  vortex  position  for  which  vortex  sink  strength  tends  to  zero.  Also, 
the  total  circulation  is  not  very  sensitive  to  chordwise  vortex  position,  but 
the  required  vortex  sink  strength  is  quite  sensitive. 

The  effect  of  moving  the  jet  chordwise  is  demonstrated  by  Figure  3«  Here  the 
jet  iB  located  at  ,40c  and  tne  vortex  at  .20c,  .25c,  and  ,30c  respectively. 
Comparing  Figures  2  and  3,  it  is  seen  that  a  chordwise  shift  of  the  jet  and 
vortex,  keeping-  the  relative  chordwiBe  spacing  about  the  same,  has  little 
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effect  on  the  overall  behavior  with  vertical  spacing.  In  every  case,  placing 
the  vortex  and  the  jet  closer  together  causes  increased  sink-strength  in  tne 
vortex,  a  result  to  be  expected  physically  since  the  jet  axial  flow  induces 
axial  core  flow  in  the  vortex. 

Figure  4  presents  a  composite  of  versus  a  data.  The  dashed  curve  is  tne 
theoretical  lift  curve  from  thin  airfoil  theory.  The  measured  for  a  flat 
plate  is  shown  as  indicated,  and  the  curve  resulting'  from  the  present  two- 
dimensional  model  is  also  3hown.  In  this  case  the  jet  is  at  .35c  chordwise 
and  .075c  vertical  and  the  vortex  is  at  ,25c  chordwise  and  at  ,125c  vertical. 
The  jet  strength  is  agam  mj  =  2.0,  slightly  lower  than  the  equivalent  value 
calculated  from  three-dimensional  experimental  data  at  =  .46.  Comparison 
of  the  theoretical  values  for  plain  potential  flow  relative  to  that  for 
vortex  augmentation  predicts,  for  example,  at  a  =  25°  a  lift  augmentation  of 
about  12$.  It  should  be  noted  that  the  lift  in  both  cases  is  computed  from 
circulation  only.  Thus,  in  the  vortex  augmented  case,  the  lift  coefficient 
is  computed  from 


(K  +  K  ). 
s  a  v' 


(6) 


Also  snown  in  the  figure  are  several  two-dimensional  values  of  lift  coefficient 
as  computed  from  experimental  results  for  =  .46.  To  arrive  at  these  values, 
the  incremental-lift  data  due  to  blowing  over  the  wing,  as  given  m  Figure  26 
of  the  experimental  part,  Volume  I  of  this  report,  was  converted  to  equivalent 
two-dimensional  values  by  use  of  the  simple  formula, 
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<g>  From  test  data,  Cj,  =  J»6  (corrected  to  2-D) 

- Potential  Theory  (Flat  Plate) 

— — Present  2-D  Theory  (Airfoil  only) 

- Present  2-D  Theory  (Vortex  +  Airfoil  Lift) 

- Helmholtz  Theory  (Ref.  12  ) 


/ 

/ 

,  / 


s 

tj 

! 


-  15  - 


These  increments  were  then  added  to  the  potential  flow  theoretical  curve 
without  spanwise  blowing.  Thus,  it  is  implicitly  assumed  that  the  vortex 
increment  for  a  flat  plate  would  be  the  same  as  for  the  three-dimensional 
airfoil  section.  The  agreement  between  theory  and  experiment  is  seen  to 
be  reasonable,  but  with  a  different  slope  in  the  data.  This  difference  in 
slope  is  probably  caused  by  the  change  in  the  actual  position  of  the  vortex 
relative  to  the  airfoil  surface  as  a  function  of  angle  of  attactc. 


It  is  observed  in  smoke  flow  results  that  the  vortex  moves  up  relative  to 
xne  airfoil  with  increasing  angle  of  attacx,  and  slightly  forward.  An  indi¬ 
cation  of  the  movement  of  the  vortex  with  angle  of  attack  can  be  obtained 
from  the  data  of  Figure  4  by  determining'  the  vertical  shift  of  vortex  position 
required  in  order  to  give  a  theoretical  value  of  lift  matching  the  values 
derived  from  test  data.  The  results  for  vertical  vortex  position  versus  a 
is  shown  in  Figure  5,  and  an  approximately  linear  relationship  with  a  for 
the  higher  angles  of  attack  can  be  seen  for  each  value  of  mj.  The 

compatible  variation  of  K  for  the  vortex  is  given  in  Figure  6.  Also, 

va 

for  the  more  optimum  jet  and  vortex  positions,  shown  in  those  figures  and 
as  indicated  by  the  test  results  of  Volume  I,  tne  vortex  sink  magnitude 
tends  to  be  on  the  order  of  twice  the  resultant  jet  source  strength. 


Jote  tnat,  for  constant  vortex  heights,  tne  magnitude  of  Kv  changes  very 


little  in  the  two-dimensional  solutions  for  a  factor  of  10  range  of  m^C^). 


This  implies  that  optimizing  the  vortex  location  is  of  great  importance 
for  minimizing  the  spanwise  blowing  energy.  The  K  *s  in  Figure  6  are 


normalized  to  a  reference  value,  K  ,  at  a  given  vortex  location, 

ref 


FIGURE  6  DERIVED  CIRCULATION  STRENGTH  CHANGES  WITH  ANGLE  OF  ATTACK 
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y  =  .125c,  and  at  an  attitude  a  =  25° »  which  is  sufficiently  hign  that 
the  flat  plate  results  should  apply  to  the  tested  wing,  since  full 
separation  would  be  expected  for  both  surfaces  without  any  blowing  momentum 
input.  Figure  5  indicates  a  further  reduction  in  yv  as  the  stall  a  is 
approached  for  a  given  airfoil  "apparent"  camber  and  thickness  increase  due 
to  the  jet  displacement  thickness  over  the  wing.  As  approaches  zero  the 
actual  incipient  stall  angle  (14  degrees)  for  the  tested  airfoil  (cf.  Vol.  I, 
Figure  29)  is  indicated,  as  expected.  For  any  value  of  the  flow  tends  to 
reattach  below  tne  natural  stall  angle,  making  discrete  vortex  formation 
unlikely  as  indicated  in  Figure  5*  There  is  also  an  upper  limit  angle  for 
each  beyond  whicn  stall  is  bound  to  occur  because  of  insufficient  blow¬ 
ing  momentum.  These  points  will  vary  with  configuration;  but  the  preliminary 
analyses  show  that  vertex  burst  at  some  station  on  the  wing  will  indicate  that 
this  augmented-’lift  stall  is  being  initiated.  Fortunately,  most  test  data 
(Refs.  3>  7  and  8)  show  that  even  past  this  point  the  maximum  lift  does  not 
decline  very  rapidly. 

Finally,  the  estimated  circulation,  K  ,  corrected  for  finite  wing  camber, 

e 

thickness  and  mT(C..)  as  well  as  incidence  angle,  can  be  obtained  by  cross 

plotting  (e.g.,  for  mj  =  2.0)  data  in  Figure  5  with  incidence, a  ,  from 

Figures  5  and  6  to  determine  the  variation  of  effective  sectional  circulation, 

Kv  ,  with  vortex  height,  y^/c  as  shown  in  Figure  7»  At  lower  nij's  the  reduced 
e 

displacement  thickness  of  the  jet  on  the  upper  surface  allows  a  smaller  forward 
airfoil  "cavity"  depth  for  a  given  angle  of  attack,  so  the  vortex  has  much 
less  strength  and  is  much  closer  to  the  surface  than  for  higher  jet  momentum 


coefficients. 


Figure  7  thus  can  be  used  to  obtain  preliminary  estimates  of  the  two- 
dimensional  Kv  and  vortex  locations  for  input  to  the  three-dimensional 
viscous  analysis  at  any  spanwise  location  in  the  vortex  feeding  sheet 
formation  region.  For  other  wing  sections  with  significantly  different 
leading  edge  radii  or  shapes,  these  predictions  at  lower  angles  of  attack 
may  not  be  very  accurate  since  the  average  height  of  the  vortex  can  be 
affected  as  discussed  in  Section  5«2*4»  Volume  I  of  this  report.  However, 
changes  such  as  in  thickness  or  camber  of  the  wing  section  may  be  approxi¬ 
mately  accounted  for  if  the  natural  (C^  =  0)  incipient  stall  angle  for  tne 
airfoil  is  known.  In  that  case,  the  assumption  is  made  that  tl  >  vortex 
height  at  high  incidence  is  not  greatly  affected;  but  tne  slope  of  the  low 
curve  in  Figure  5  changes  to  conform  with  the  intercept  on  the  abscissa, 
which  should  be  the  stall  angle  for  the  new  airfoil. 
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3.0  THREE-DIMENSIONAL  THEORY  FOR  VORTEX-LIFT  AUGMENTATION 

3.1.0  INTRODUCTIONJPO  THE  THRE&-D IMEN S ION. AL  ANALYSIS 

The  analysis  of  cpanwise  blowing  for  vortex-lift  augmentation  must  ultimately 
resolve  the  problem  of  regions  of  high  energy  jet  and  vortex  flows  in  close 
proximity  to  a  wing  surface,  as  characteristic  of  many  powered  lift  systems. 
Ideally,  a  mathematical  model  should  include  three-dimensional  boundary 
layer  interactions  with  the  jet  and  vortex  flow  field.  However,  the  present 
research  objective  was  to  obtain  a  physical  understanding  of  the  various 
flow  phenomena  involved.  Detailed  pressure  measurements,  flow  visualizations 
of  the  leading  edge  and  trailing  vortex  systems  and  LDV  measurements  of  a 
jet  issuing  tangent  to  a  plate,  were  used  to  gain  an  insight  into  the  flow 
mecnanisms  required;  but  at  present  the  relative  importance  of  eacn  of  these 
flow  fields  can  be  established  only  by  analysis  until  more  detailed  measure¬ 
ments  are  made,  particularly  concerning  mixing  in  a  cross  flow.  In  the 
analytical  model,  the  wing,  the  augmented  leading  edge  vortex  with  a  feeding 
sheet,  and  the  jet  centerline  are  represented  by  singularities  consisting  of 
source-sinks  and  vortex  (lattice)  distributions.  Appropriate  boundary 
conditions  are  imposed  to  determine  the  singularity  strengths  on  the  surface. 

The  principal  empiricisms  required  are  associated  with  the  trajectory  and 
enhancement  of  the  spanwise  blowing  jet  imbedded  in  the  decelerated  cross¬ 
flow  region  over  the  leading-edge  of  the  wing  upper  surface.  These  formulas 
are  taken  with  minor  modifications  from  existing  literature  and  discussed  in 
Appendix  A. 
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3.2.0  aYTOTHE3ISJ^R_STi£ILIZATION__Of  THE  LEAI)HfG-EDGE_VORTEX 

The  lift  augmentation  phenomena  described  herein  are  most  pertinent  to  wing 
angles-of-attack  where  highly  separated  flow  from  the  leading-edge  is  normally 
present.  Experimental  evidence  is  sufficient  to  prove  that,  with  spanwise 
blowing,  a  steady  vortex  will  form  in  the  lee  of  the  separation  region;  so 
tne  analysis  can  proceed  without  further  comment.  However,  an  hypothesis  can 
be  made  for  the  mechanism  of  leading-edge  vortex  stabilization,  based  on  purely 
analytical  arguments.  Without  rigorous  proof,  a  brief  argument  is  presented 
in  the  interest  of  possible  better  understanding  of  vortex  stability. 

Without  spanwise  blowing,  at  higher  Reynolds  Numbers,  a  wing  experiences  buffeting 
caused  oy  quasi-periodic  shedding  of  vorticity  and  breakdown  into  turbulence 
in  the  near  wake.  The  methods  of  Refs.  1  and  6  have  been  proven  to  be  capable 
of  predicting  the  dominant  spectral  frequencies  in  the  wake.  In  particular, 
the  Strouhal  frequency  and  its  harmonics  are  derived  in  closed-form  using  tne 
Uavier-Stokes  equation  as  summarized  from  Reference  1  in  Appenuix  B.  The  solution 
is  shown  to  be  largely  dependent  upon  the  Rossby  number,  Rq,  (ratio  of  linear- 
to-angular  momentum)  and  the  Radial-inflow  Reynolds  number,  R^.  Prom  the 
form  of  the  equation  of  fluid-dynamic  motion,  it  is  apparent  that  stable,  but 
periodic  vortex  modes  are  predicted  when  R^  is  negative-  so  that  the  outer 
vortxcal  layers  of  the  separation  region  supply  net  radial  inflow  in  a  sustained 
manner.  Under  natural  conditions,  the  vortex  Will  continue  to  grow  until  the 
core  axial  flow  accelerates  to  a  burst  condition  (Ref.  13),  or  until  the 

wake  pressure  defect  can  no  longer  restrain  the  vortex,  so  that  a  fully 
separated  wake  forms  with  periodically  shed  vorticity. 


-  23  - 


However,  spanwise  blowing  is  known  to  prevent  these  previously  described 
oscillatory  conditions.  Tc  understand  this,  one  must  recognize  that  the 
transformed  equation  (Appendix  B)  will  indicate  stable  formation  and  axial 
growth  patterns  if  R^.  is  positive.  This  is  only  possible  in  the  upper 
part  of  tne  "cavity"  flow  region,  as  sketched  in  Figure  8,  where  the  outer 
vortical  layers  of  the  leading-edge  feeding  sneet  begin  to  swirl  in  towards 
the  vortex  core.  However,  simultaneously  the  co-rotating  jet  induces 
velocities  which  oppose  those  just  behind  the  vortex.  This  reverses  the 
entrainment  process  at  the  vortex  periphery,  and  could  be  interpreted  as 
vortex  radial  outflow.  The  mechanism  for  stabilizing  the  outer  vortical 
layers,  and  preventing  the  periodic  form  of  vortex  shedding  is  thereby 
provided.  As  discussed  in  the  two-dimensional  analysis,  the  jet  also  acts 
as  a  source,  providing  flow  toward  the  inner  core  of  the  vortex.  This  two¬ 
fold  character  of  the  jet  is  also  described  in  Section  Volume  I  by 

using  test  observations  on  the  high  spreading  rate  at  the  core  of  the 
spanwise  jet.  These  lateral-spread,  wall- jet  velocities  combine  with  the 
feeding  sheet  elements  over  the  wing  surface,  as  observed  in  Figure  9(a), 
to  entrain  mass  and  momentum  into  the  underside  of  the  vortex.  The  jet 
thus  augments  the  vortex  core  strength  in  a  two-dimensional  sense.  In  tests, 
however,  it  has  been  observed  that  excessive  jet  velocity  can  blow  the  vortex 
off  the  wing.  Therefore,  an  equilibrium  must  be  established  wherein  the 
detrainment  (outflow)  from  outer  vertical  layers  must  balance  the  entrain¬ 
ment  (sink  effect)  at  the  vortex  core.  Three-dimensionally,  then,  it  is 
probable  that  there  is  a  net  reduction  in  the  amount  of  fluid  which  would 
otherwise  accelerate  ir  to  the  core  and  spiral  out  the  vortex  axis.  This  also 
explains  the  delay  in  vortex  burst,  caused  by  spanwise  blowing,  until  the 
vortex  curves  aft  to  coalesce  with  the  flap  outer-edge  or  wing-tip  vortex. 


OBSERVED  PHYSICAL  DESCRIPTION  OF  VORTEX  LIFT  CONTROL 
(With  full  leading  edge  separation  ) 


VORTEX  LIFT  CONTROL  FROM  TEST  FLOW  VISUALIZATION 
AT  a  *  1*1°  (REF.  3) 
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The  veriation  of  vortex  strength  caused  by  increased  entrainment  in  the 
mixing  region  is  tractable  by  a  method  similar  to  that  of  Owen  (Ref.  14) » 
where  the  hard  core  is  bounded  by  a  finite  annulus  characterized  by  a  time- 
dependent  turbulent  shear  stress  and  the  cascading  of  energy  with  frequency. 
This  will  be  discussed  in  more  detail  in  a  later  section  on  three-dimensional 
stability  and  viscous  effects. 

The  form  of  equation  used  in  Reference  1  for  axial  motion  is  probably  not 
the  most  suitable  for  aperiodic  flows  where  the  leading-edge  vortex  is 
growing  rapidly  and  approaching  burst  conditions.  Also,  the  axis  of  the  jet, 
and  hence  the  vortex  core,  assumes  a  rearward  trajectory  on  the  wing  surface 
because  of  the  high  local  crossflow  velocity.  The  general  problem  of  a  free 
jet  penetrating  a  cross-flowing  stream  has  received  much  study  recently 
(Refs.  15,  16  and  17),  but  modifications  to  account  for  skewed  impingement 
tangentially  on  a  lifting  surface  are  derived  semi-empirically  herein  since 
no  suitable  analytic  methods  appear  to  be  available. 

3.5.0  ORSISVAT ION S_0N  THE  ^-Y0RTEX_roTER4CT  IQja 

3.5*"1  background  Analyses  and  Tests 

Reference  to  Figure  8  which  was  deduced  from  current  experiments,  shows 
in  more  detail  how  the  primary  vortex  and  the  jet  interact  beneficially. 
Although  some  approximate  solutions  are  available  for  either  a  jet  in  a 
cross-flow  (e.g.,  Refs.  15  and  18)  or  for  a  viscous  vortex  forming  in  a  near¬ 
wake  (Refs.  1  and  5)  or  leading-edge  separation  region  (Ref.  6),  no  solutions 
are  available  for  a  combination  of  these  two  issuing  coaxially  from  a  wall 
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(fuselage)  and  parallel  to  a  surface  (wing  or  flap).  Tne  Laser  Doppler 
Velocimeter  (LDV)  measurements  shown  in  Figure  11  of  Volume  I  indicate  that 
the  jet  velocity  profile  behaves  much  like  a  wall  jet  past  the  under-expanaed 
region  where  Coanda  type  reattacnment  occurs.  Visual  evidence  consisting  of 
photographs  of  the  jet  erosion  pattern  over  the  upper  surface  during  the  LDV 
tests,  reveais  that  the  jet  spread  angle  Y  is  about  4  times  larger  tnan  that 
for  a  free  jet,  largely  due  to  increased  entrainment  due  to  mixing  with  tne 
rather  thick,  decelerated  inner  region  of  the  boundary  layer.  A  purely 
theoretical  approach  to  describe  such  a  complex,  mixed  flow  region  is  beyond 
the  scope  of  the  present  work.  Detailed  numerical  solutions  for  wall-jets  may 
De  conceptually  possible  but  the  cross-flow  and  vortex  and  bounaary-layer 
viscous  terms  would  require  a  higher  order  theory.  For  the  present  purpose, 
a  more  intuitive  formulation  of  the  problem  is  employed  to  test  the  validity 
of  the  hypothetical  flow  mecnanisms  above  the  wing  surface.  Direct  measure¬ 
ments  to  substantiate  the  velocities  in  this  region  must  await  tests  using  a 
backseat tering  laser-dopp'ler  or  similar,  non- interfering  measuring  system. 

Flow  visualization  has  been  helpful,  but  only  the  surface  pressure  measure¬ 
ments  are  quantitative.  Thus,  a  semiempincal  descriptive  analytical  model 
is  used,  calling  upon  available  formulae  whenever  possible.  The  flow 
augmentation  observed,  however,  indicates  that  the  coefficients  in  such 
formulae  are  likely  to  need  moaif ication. 

A  literature  survey  has  been  made  on  the  aforementioned  topics,  and  a  summary 
of  those  semi-empirical  formulae  applied  herein  is  given  in  Appendix  (a). 
Detailed  justifications  for  each  equation  selection  are  not  always  given 
since  in  many  cases  alternative  expressions  are  available,  but  the 


differences  in  numerical  results  from  a  variety  of  references  are 
considered  minor,  in  most  cases. 

3.3.2  Jet  Trajectory 

Two  primary  assumptions  are  made.  First,  the  presence  of  an  image-plane 
due  to  the  jet  over  the  wing  upper  surface  would  cause  some  jet/vortex 
cross-sectional  distortion  (e.g. ,  Ref.  18),  but  this  is  considered  to  be  a 
higher  order  effect,  so  the  variation  of  maximum  axial  local  velocity,  Uffl, 
and  the  trajectories  of  a  circular  jet  issuing  from  a  vertical  wall  in  a 
cross  flow  are  assumed  to  apply  in  basic  form.  It  is  recognized  that  sucn 
a  jet  forms  a  plume  consisting  primarily  of  a  pair  of  counter-rotating 
vortices  which  curve  towards  the  mainstream  direction.  However,  it  is  observed 
in  the  case  of  spanwise  blowing,  conceptually  shown  in  Figures  8  and  b,  that 
only  one  swirling  jet  issues  in  close  proximity  to  the  wing  upper  surface, 
but  curving  aft  behind  the  leadnug  edge  vortex.  As  the  jet  profile  develops 
in  the  spanwise  direction,  its  characteristics,  revealed  by  the  preliminary 
LDV  measurements,  approach  those  of  a  wall- jet.  The  upstream,  spreading 
component  of  velocity  is  terminated  by  a  dividing  streamline,  separating  the 
jet  and  vortex,  yet  revealing  a  path  for  the  jet  lateral  velocity  components 
to  provide  mass  flow  into  the  core  of  the  primary  vortex. 

Thus,  for  a  reasonable  fluid  dynamic  analogy,  the  assumption  must  be  made 
tnat  the  counter-rotating  partner  of  the  swirling  jet  is  imbedded  in  an 
image  plane  formed  by  the  wing  upper  surface  as  sketched  in  a  cross-section 
of  Figure  8. 
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3.3*3  The  Leading-Edge  Vortex 

For  a  finite  wing  leading  edge  at  high  angle  of  attack,  a  feeding  sheet 
forms  somewhat  as  sketched  in  Figures  8  and  9(a).  The  data  of  Reference  3 
gives  substantiation  of  the  assumption  that  the  vorticity  emanating  from 
the  wall  (fuselage)  region  of  jet  formation  curves  forward  over  the  upper 
surface,  under  the  jet  and  "unbound”  vortex  and  then  is  abruptly  turned 
upwards  at  a  fairly  sharp  angle.  The  curved  sheet  is  approximated  by  a 
straight  vortex-element  sheet  as  sketched  in  Figure  9(b).  The  vortex  elements 
are  assumed  to  pass  from  the  leading-edge  to  the  primary  vortex  center  m  tne 
manner  assumed  by  Brown  and  Michael  (Ref.  19).  However,  much  improvement 
in  results  is  obtained,  e.g. ,  for  delta  wings,  using  this  greater  detail  of 
the  formation  region  over  the  wing  surface.  A  sample  comparison  of  this 
method,  for  *  0,  with  that  of  Polhamus  (Refs.  20  and  21)  for  an  aspect 
ratio  2.0  delta  wing  is  shown  in  Figure  10.  The  agreement  in  lift  partially 
substantiates  the  selected  feeding  sheet  arrangement. 

For  a  surface  of  infinite  extent,  a  wall-jet  tends  to  accelerate  in  the 
favorable  pressure  gradient  created  under  a  vortex  of  proper  sign,  as  exists 
in  Figure  8,  but  it  may  not  greatly  affect  the  strength  of  the  vortex  above 
the  surface.  The  test  data  (Figure  10  of  Volume  I)  indicate  that  normally 
the  jet  spreads  at  a  much  greater  rate  than  does  a  free  jet,  probably  due  to 
surface  friction  effects  and  entrainment  with  a  highly  turbulent  boundary 
layer.  In  this  case,  however,  the  presence  of  a  strong  feeding  sheet  at  the 
wing  leading  edge  will  decelerate  the  forward  lateral-spread  component  of 
the  "wall-jet"  abruptly,  eventually  resulting  in  a  streamwise  separation 
profile  near  the  surface  and  a  conversion  of  its  vorticity  in  the  process. 

Some  acceleration  of  the  feeding  sheet  velocity  must  occur  as  the  jet  velocity 
is  turned  upward  approximately  normal  to  the  free  stream  direction. 


FIG.  9(b)  LEADING-EDGE  VORTEX  LIFT -CONTROL  WITH  VORTEX- LATT I  CE  AERODYNAMICS 
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FIG.  10  COMPARISON  OF  LEADING-EDGE  VORTEX  CONTRIBUTION  TO  TOTAL  LIFT  USING  A  DELTA  WING 
CASE  FROM  PRESENT  THEORY  FOR  COMPARISON  WITH  THAT  FROM  POLHAMUS'  THEORY  FOR 
VARIOUS  FEEDING  SHEET  AND  POTENTIAL  THEORY  ASSUMPTIONS  ON  NORMALIZATION. 
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Concurrently,  as  described  in  Polhamus'  leading-edge  suction  analogy,  Refs. 
20  and  21,  the  thrust  force  of  this  turning  across  the  feeding  sheet  region 
results  in  momentum  exchange  causing  a  lift  increment  of  equal  magnitude. 

The  vorticity  involved  is  converted  to  added  circulatory  lift  by  augmenting 
the  strength  of  the  leading-edge  vortex  suspended  in  the  low  pressure  region 
ahead  of  the  jet  and  behind* the  feeding  sheet.  The  simplifying  assumption 
is  made  that  only  the  positive  (lift-contributing)  increment  of  the  lateral 
jet  profile  circulation  need  be  considered  in  estimating  the  vortex 
augmentation,  AFj. 

3-3*4  Secondary  Vortices 

Presumably  the  negative  vorticity  components  contribute  to  the  formation 
of  many  turbulent  secondary  vortices  in  a  highly  viscous  region  very  near 
the  base  of  the  feeding  sheet.  Here,  wind  tunnel  tests  indicate  likely 
presence  of  one  or  more  small,  counter-rotating  secondary  vortices  forward 
of  the  primary  lifting  vortex  as  sketched  in  the  Bectional  view  of  Figure  8. 
Probably,  a  series  of  such  minute  eddies  dissipates  the  energy  of  the  lower 
side  of  the  wall  jet  profile. 

3.4*0  W AIJrvJET JOM  PROPAGATION  AND  DKAY 

3*4*1  Leading-Edge  Vorticity  Estimation  from  Jet  Profile 
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The  LDV  measurements,  Figure  13  of  Volume  I,  indicate  that  the  part  of  the 
wall-jet  profile  Um(x,y, z)  above  the  maximum  velocity  point,  Ayffl,  (adapting 
terminalogy  from  Ref.  17)  can  be  approximated  by  a  straight  line.  Using 
the  coordinate  system  in  Figure  A.l  of  Appendix  A,  commonly  accepted  for  a 
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jet  xn  a  cross-flowing  stream  (Refs.  15  and  16),  the  linear  portion  of 
the  lateral  jet  velocity  profile  just  beneath  the  vortex  axis  may  be 
expressed  approximately  as 


dw 

dy 


tan  7  cos  0  . 

i 


(8) 


v  . 
ni 

where:  (Refer  go  sectional  sketch  in  Figure  9(b)) 

zv  designates  the  chordwise  location  of  the  nta  link  of  the 
ni 

primary  vortex  at  the  ith  spanwise  section  (at  x.). 


Uy  is  the  jet  axial  velocity  component  at  the  above  vortex  location 
ni 

(see  Appendix  A). 


7  is  the  jet  spread  half-angle  (Figure  9(b)). 


9^  is  the  trajectory  angle  of  the  jet  centerline,  relative  to  the 
cross-stream  direction,  at  section  i. 
and 

is  the  effective  depth  of  jet  beneath  the  vortex  at  the  dividing 
streamline  (see  Figures  8  and  9(b)),  above  the  level  Aym> 


Since  no  data  are  available  for  a  jet  with  these  boundary  conditions,  the 
LDV  measurements  (Figure  11  of  Volume  I)  are  used  to  justify  a  value  of 
2  1.7  Do  or  approximately  a  constant  times  the  nozzle  diameter,  Dq,  as 
sufficiently  accurate  since  the  unaccelerated  wall  jet  profile  is  vanishing 
near  Ay  =  <5^.  Note  that  Dq  is  an  effective  diameter  (Appendix  A). 

Finally  the  vorticity  must  be  integrated  over  each  interval,  Az  , 

n* 

from  the  wing  leading-edge  to  the  vortex  chordwise  position. 
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Az  =  z  -  x.  tan  A 
v  .  v  .  1 

m  ni 


to  yield  the  wall-jet  vorticity  available  for  vortex  augmentation. 


3«4«2.0  Vortex  Augmentation  and  Decs 


3 . 4 • 2 . 1  Augmentation 


With  q*  as  the  vector  velocity  across  a  line  element  ds  the  contour  integral 
for  the  jet  augmentation  of  vortex  circulation,  AI^.  due  to  viscous  spreading 
in  the  boundary  layer  at  a  section  transverse  to  the  vortex  centerline  may 
be  expressed  as  a  Stokes’  integral: 

AT  =  ffi  curl  q*  •  ds 


/  7  ( 


bv  _  bv 
by  ~  bz 


Within  the  accuracy  of  the  jet  profile  data  on  the  chordwise  gradient  in 
proximity  to  the  stagnation  streamline  (Figure  8)  the  assumption  may  be 
made  that  bv  ^  bw  .  For  the  higher  blowing  momentum  coefficients,  the 

5z  by 

bounding  and  stagnation  streamlines  remain  at  a  nearly  constant  position 
relative  to  the  wing  leading  edge.  Thus,  the  jet  is  restrained  from 
penetrating  any  farther  forward,  and  it  must  expend  its  energy  in  entrain¬ 
ment  with  the  vortex  above  the  surface.  The  dominant  gradient  for  this 

effect  is  bv  although  civ  has  implications  with  respect  to  vortex 

by  £z 

stability  (burst)  as  will  be  discussed  later. 

Substitution  of  equation  (8)  into  (10),  with  dv  ~  o,  yields 


Within  the  accuracy  required  for  this  purpose,  it  is  assumed  that  Uy 

ni 

is  approximately  constant  over  each  chordwise  increment  of  integration, 

Azni,  principally  influenced  by  only  the  nto  primary  vortex  link.  Thus,  to 

simplify  the  integration,  U  and  Q ^  are  not  considered  as  functions  of 

ni 

y  and  z,  particularly  for  the  cases  of  most  significance  where  the  primary 
vortex  is  formed  very  near  the  leading-edge  feeding  sheet.  The  result  may  be 
expressed  in  the  form, 

ATj  s  Uy  (zv  -  x^anA)  tanTcos^)  (12) 

ni  ni  ni 


This  equation  provides  an  estimate  of  the  increase  in  circulation  above 
that  predicted  by  the  2-dimensional  theory  (2ttK)  which  is  used  as  the 
starting  value  for  the  vortex  at  the  first  inboard  wing  station  beyond  the 
under-expanded  jet  distance,  x  ,  from  the  fuselage  sidewall  or  nozzle.  At 
each  spanwise  station,  x^, the  cumulation  of  circulation  increments  from  the 
primary  vortex  links,  given  by  equation  (12)  is  added  to  the  two-dimensional 
values. 


II 


5 

!  s 


3 . 4 . 2 . 2  VortjsxJDecay 

It  was  found  that  for  the  above  rather  high  vortex  strengths  the  vortex 
"ages"  rather  rapidly  according  to  Professor  Owen’s  theory  (Ref.  14),  as 
described  in  Appendix  A.  Previous  successful  experience  with  this  metnod  over 
a  wide  Reynolds  Number  range,  (Ref.  25),  encouraged  its  application  herein. 

The  time  varying  turbulent  eddy  viscosity  resulting  therefrom  decays 


-  36  - 


the  vortex  strength  in  a  somewhat  linear  manner  to  values  typically  1/3 
of  those  predicted  by  the  two-dimensional  potential  theory.  An  iterative 
calculation  procedure  is  thus  indicated;  but  it  was  found  that  the  rule-of- 
thumb  first  estimate  of  the  reduced  vortex  strength  and  trajectory  vari¬ 
ations  resulted  in  convergence  for  practical  purposes  in  one  or  two 
iterations. 

30.0  THE  VORTK-UTTICE/SOUEC^SIliK_PROGRAl'I 
3.5.1  The  Computer  Simulation  Model 

The  above  equations,  along  with  the  semi-empirical  formulations  for  jet  ana 
vortex  growtn,  entrainment,  and  viscous  dissipation  found  in  Appendix  A, 
provide  the  basic  elements  necessary  to  simulate  vortex-lift  control  as 
depicted  in  Figure  b(b).  However,  the  application  of  such  elements  in  a 
vortex  lattice  program  (e.g.  Refs.  15  and  22-24)  for  powered/augmented  lift 
analysis,  as  delineated  in  Figure  11,  is  more  complex  and  extensive  than  was 
originally  intended  in  this  research  program.  The  phenomena  are  so  highly 
three-dimensional  and  viscous  that  the  limitations  of  the  two-dimensional 
theory  became  evident  early  in  the  analysis.  Nevertheless,  the  2-3)  theory 
is  essential  to  provide  starting  conditions  for  logical  vortex  locations 
and  strengths  and  source/sink  effects. 

Ideally,  the  singularity  distributions  from  the  jet,  vortex  feeding  sheet 
and  lifting  surface  should  all  be  included  simultaneously  in  the  computation 
of  mutual  induction  effects  while  satisfying  the  flow  tangency  boundary  con¬ 
ditions  at  collocation  points  throughout  the  wing.  This  objective  is 
attempted  in  the  final  steps  as  shown  in  Figure  11;  but  certain  non-linearities 
present  in  such  a  viscid/ mviscid  interaction  process,  particularly  in  the 
event  of  vortex  burst,  imply  tiiat  assumptions  of  superposition  must  be  made 


INITIAL  SINK 
ESTIMATE  FROM 
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FIGURE  11.  VORTEX-LIFT  CONTROL  FLOW  DIAGRAM  OF  THREE-DIMENSIONAL 

VISCID  AND  I NV I  SC  ID  COMPUTER  PROGRAMS 
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at  intermediate  steps  in  the  program  calculations.  Hopefully,  iterations 
on  all  or  part  of  the  steps  in  Figure  11  compensate  and  lead  to  minimum 
errors.  For  example,  the  jet  trajectory  is  dependent  upon  the  mean  upper 
surface  crossflow  velocity.  This  is,  in  turn,  related  to  the  vortex  strength 
variation  since  the  circumferential  velocity  creates  a  backflow  component 
tangentially,  close  to  the  surface  (e.g. ,  see  the  Figure  8  typical  cross-section), 
rut  tne  vortex  strength  depends  on  the  jet  strength  and  entrainment  variation 
in  the  boundary-layer/,,wall-jet"  region;  so  an  iteration  is  indicated  in  the 
flow  diagram  of  Figure  11  with  a  convergence  check  on  the  jet  and  vortex  paths. 

If  tne  leading-edge  vortex  strengths  vary  significantly  from  those  for 
prevxous  steps,  the  feeding  sheet  is  modified  accordingly,  and  new  induction 
effects  are  calculated  to  enforce  the  boundary  conditions  again. 

Fortunately,  the  pressure  coefficients  did  not  vary  significantly  on  the 
average,  so  few  iterations  were  required  for  the  optimum  cases  studied. 

In  this  respect,  the  program  is  not  sensitive  to  the  distributed  sink  strengths 
except  at  low  attitudes  wnere  close  proximity  of  the  unbound  vortex  elements 
to  the  wing  surface  induces  some  locally  high  velocities.  Such  cases  are 
not  optimum  for  vortex  lift  augmentation,  and  thus  should  not  be  of  great 
concern.  However,  the  problem  waB  found  to  be  much  alleviated  by  applying 
smoothing  techniques  or  spline  curve-fitting  techniques  to  the  most  severely 
fluctuating  induced  velocities  in  the  previous  iteration  before  proceeding 
to  the  next  simultaneous  equation  solution.  Such  methods  are  recommended 
(Refs.  27  and  28)  where  finite-element,  singularity-distribution  methods 
are  employed.  The  solution  convergence  is  much  more  rapid  for  optimum 
smoothing  (Ref.  29),  and  the  mathematical  possibility  of  introducing 
non-unique  solutions  is  minimized. 
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In  the  future,  a  solution  is  conceptually  possible  by  a  more  direct 
iterative  approach  combining  the  time/spacially  varing  lattice  of  Ref.  15 
with  t jne-dependent  turbulent  boundary  layer  methods  (Ref.  30);  but  the  semi- 
empirical  method  proposed  herein  is  more  feasible  at  this  time  since  it 
employs  approximate,  closed-form  expressions  for  viscous  and  entrainment 
processes  along  with  the  well  known  and  proven  vortex-lattice  tecaniques. 

More  detailed  descriptions  of  the  latter  programs  are  given  in  References 
15,  and  22  through  24. 

3-5.2  Input  Conditions  for  the  3-D  Solution 

The  program  sequence  is  shown  in  the  flow  diagram  of  Figure  11.  The 
coordinate  locations  for  each  panel  of  the  wing,  flaps,  or  other  portions  of 
the  airplane  are  computed  automatically  in  the  vortex-lattice  (V-L)  computer 
program  by  merely  specifying  the  corner  points  of  each  major  section  which 
does  not  have  either  chordwise  or  spanv>°s  discontinuities,  Thu6,  it  is 
convenient  to  first  make  a  V-L  run  at  the  desired  attitude, o,  to  obtain  the 
coordinates  of  the  load  points,  through  which  the  spanwise  bound  vortex  link 
of  each  panel  passes,  and  the  collocation  (control)  points  where  the  exact 
tangency  boundary  conditions  are  imposed  in  the  matrix  solution.  This 
provides  a  check  against  the  unpowered  lift,  or  potential  flow  solutions. 
Also,  these  coordinates  are  required  as  input  for  the  solution  of  tne  viscous 
jet-vortex  (J-V)  equations  (Appendix  A)  in  a  cross  flowing  stream. 

The  chordwise  coordinates  and  initial  strengtns  of  the  vorte-:  and  jet  for  a 
given  C ,  and  a  axe  specified  from  the  two-dimensional  solution  at  a  char¬ 
acteristic  spanwise  station.  This  may  be  either  at  the  wing  root  or  at  any 
position  for  outboard  spanwise  blowing.  The  jet  and  vortex  velocities  and 
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trajectories  are  computed  for  as  many  nozzles  as  are  present.  The  jet 
augmentation  (equation  12)  and  viscous  dissipation  from  Owen's  theory, 

Ref.  14,  are  accounted  for  in  obtaining  an  improved  vortex  strength  variation, 


Kni,  witii  span  which  is  tnen  input  to  the  V-L  program,  properly  accounting 


for  the  cross-flow  trajectory  of  the  vortex  centerline.  As  discussed  previous¬ 
ly,  the  unbound,  locked  vortex  feeding  sheet  is  assumed  to  emanate  from  the 
fuselage  wall  and  boundary-layer  region  by  means  of  equally-spaced  discrete 
vortex  elements  which  are  assumed  to  be  bound  in  the  surface  parallel  to  the 
local  sweep  angle,  and  coincident  with  the  nearest  quarter-panei-chord  bound- 
vortex  leg  to  minimize  potential  interactions  of  two  singularities  in  proximity 
to  each  other  (e.g. ,  see  Ref.  15).  These  vortex  links,  having  strengths 
consistent  with  the  previously  determined  spanwiBe  F  variation,  are  then 
assumed  to  extend  spanwise  until  they  are  turned  forward  by  the  backflow 
velocity  components  approximately  beneath  the  jet  or  locked-vortex  centerlines, 
as  appropriate  from  the  wing  geometry  and  J-V  trajectory  printout.  They  are 
then  swept  forward  to  the  wing  leading  edge  and  thence  back  up  to  the  unbound 
vortex  to  simulate  the  Leading  edge  feeding  sheet  represented  somewhat  in  the 
manner  of  Brown  and  Micheal,  Ref.  19.  Continuity  of  vortex  lines  is  ensured 
by  allowing  some  of  the  local  vorticity  to  be  shed  into  the  wake  at  an  angle 
01  /2,  based  on  experience  with  V-L  methods,  Ref.  22-24,  and  the  remainder  to 
propagate  out  the  wing  span  and  aft  to  eventual1 y  merge  with  the  wing  tip 
vortex. 


The  additional  entrainment  expected  from  the  finite  wing  leading  edge 
conditions  seems  to  be  derived  naturally  from  the  two-dimensional  solution  for 
the  range  of  practical  attitudes  and  momentum  coefficients  shown  in  i’igure  5. 


■m 


i 
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On  the  average,  the  vortex  sink  strength  varied  less  than  plus  or  minus 
25  percent  from  the  relationship 

«■  -  2  mj  (13) 

and  the  variation  above  stall  was  not  consistently  increasing  or  decreasing 
with  a  for  tne  range  of  momentum  coefficients,  .046^0^.46.  Thus,  within 
the  accuracy  of  the  two-dimensional  assumptions,  equation  (13)  was  used  as 
sufficient  for  starting  conditions  in  the  V-L  program.  The  J-V  program  computes 
values  of  d^dx,  volumetric  mass  flux  change  with  wing  span  as  discussed  in 
Appendix  A,  and  these  can  be  used  to  obtain  more  accurate  values  for  mj  if 
this  is  desired. 

n'or  a  given  ,  values  of  Kv  versus  a  can  be  established  for  constant  jet  and 

vortex  locations  which  are  very  much  like  those  shown  in  Figure  6.  However, 

based  on  earlier  discussions  in  Figure  6  it  is  also  necessary  to  account  for 

the  reduction  in  Ky  with  lesser  values  of  yy  as  shown  in  Figures  2,  3,  and  7. 

as  angle  of  attack  decreases,  the  vortex  clings  closer  to  the  upper  surface, 

and  the  separation  (cavity)  region  is  smaller,  so  vortex  strength  reduces 

typically  as  shown  in  Figure  7»  The  usual  procedure  is  v.  select  an  attitude 

sucn  as  25  degrees,  sufficiently  above  natural  stall  for  a  given  airfoil 

that  the  vortex  strength  (then  called  K  )  and  height  determined  from  the 

Ref. 

flat-plate  theory  will  be  adequate  for  a  fully  separated  flew-.  Then  tne  two 
factors  from  curves  such  as  those  in  Figures  6  and  7  are  both  applied  to 
unnormaiize  the  Kv  and  yv  for  lower  values  of  a  . 

The  3-b  program  has  the  option  01  computing  tne  jet  velocities  out  of  a 
singularity  distribution,  but  the  "lumped  parameter"  distribution  of  these 
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sources  was  not  available  in  an  economical  manner  from  any  theory  or  tests; 
so  an  alternative  procedure  was  devised  using  jet  vector  "supervelocitj.es" 
from  Appendix  A  formulas,  multiplied  by  an  empirical  factor  prior  to  application 
at  grid  collocation  points,  as  discussed  in  the  next  section. 


3o»3»1  ^uMnt.ry_  £i_l£ad_s_C£iirDut_a_t i£n_I'Ie thod^  -  ‘The  purpose  of  the  tnree- 
dimensional  program  is  to  combine  the  effects  of  jet-vortex  interactions, 
along  with  vortex  persistence  and  stability,  to  obtain  the  resultant  velocities 


over  an  arbitrary  wing  surface  represented  by  a  vortex  lattice.  Tnus  it  is 


possible  to  compute  the  incremental  pressure  coefficients,  Ac^,  associated  witn 
the  jet  and  vortex  source/sink  systems,  and  with  the  vortex  lattices  for  the 


wing  and  the  feeding  sheet. 


In  general,  tne  total  force  at  a  local  station  on  the  wing  according  to 
potential  theory  can  be  expressed  in  the  form, 

Total  Lift  =  p  (vy_L  +  vg)  (  rv_L  +  ry)AL  (14) 

where, 

V?  ^  are  velocities  which  lie  in  vertical  streamwise  planes  and 

are  induced,  at  the  load  points,  by  the  basic  vortex  lattice  systia 
system  in  the  wing. 


Vc  are  superposed  velocities,  in  streamwise  vertical  planes, 
caused  by  feeding-sheet  elements,  by  sources  and  by  sinks 
in  the  field.  (There  is  also  a  fuselage-body  thickness 
correction  included  (Ref.  31 )»  in  most  cases.) 


ry_L  are  the  local  circulations  determined  after  all  mutually- 
induced  and  all  supervelocities  have  been  applied  at  the 
collocation  points.  The  effects  of  the  feeding  sheet 
are  included  explicitly. 

F  are  the  local  circulations  in  the  wing  which  join  the 
leading-edge  feeding  sheet  and  ultimately  form  the 
separated  vortex  above  the  wing.  These  eventually  trail 
downstream  to  infinity. 

and, 

AL  is  the  cross-stream  projected  length  of  a  given  bound 
vortex  element. 

A  difficulty  is  created  by  the  cross-product  terms  of  equation  (14),  since 
these  imply  tnat  simultaneous  solution  of  the  inviscid  (V-L)  and  viscid  (J-V) 
equations  is  required  to  obtain  mathematically  consistent  loads.  The 
assumption  is  made  that  the  iterative  convergence  tests  are  sufficient  to 
satisfy  the  necessary  compatibility  and  uniqueness  conditions. 

50*3.2  £emi_2empi^ri_cal_Modi_fi_catiions_t£  Surface  £u£ery£l_oc_it_ie_s  for  Reduced 
Crossflow  -  Experience  with  vortex-lattice  methods  reveals  that  the  matrix 
of  influence  coefficients  is  adequately  determined  for  most  conventional  wings 
by  consideration  of  only  the  potential  velocity  field  beyond  the  boundary 
layer.  At  most,  a  classical  thickness  and/or  camber  correction  may  be  necessary 
to  account  for  an  effective  wing  thickness  increase  due  to  boundary-layer 
displacement  thickness.  In  the  case  of  "boundary  layers"  consisting  of 
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thick  wall- jet  profiles,  where  velocity  defects  are  significant  at  large 
distances  from  the  wing  surface,  available  literature  doeB  not  indicate  any 
simple  correction  that  has  been  successfully  tried.  This  situation  is 
aggravated  by  the  tangential  velocity  gradients,  created  by  the  unbound 
"leading-edge"  vortex,  vrn..uu  marge  with  the  jet  velocities  as  shown  in 
Figures  8  and  9  at  even  greater  distances  in  the  separated  flow  region. 


If  an  effective,  mean  upper-surface  (cavity)  flow  velocity  is  defined  as 
V  ,  then  we  may  express  the  cross-flow  velocity  ratio  for  the  jet  as 

H  -  Ve  (15) 

It  already  has  been  shown  (Equation  5  and  Refs.  6,  10  and  11)  that  the 
velocity  in  this  region  may  be  as  little  as  10  percent  of  tne  free  stream 
velocity,  V^  .  Thus  for  a  sonic  jet,  as  found  to  be  optimum  in  Volume  I 
of  this  report, 

R  »  1  (16) 

A  similar’  argument  can  be  made  for  the  ratio  of  the  effective  "streaming- 
flow"  component,  V  ,  divided  into  the  maximum  wall-jet  velocity  at  station  i 
in  the  laterally  spreading  direction  just  beneath  the  n^  unbound  vortex 
element,  given  as 


*VniC0S^i  't:an^ 


(17) 


In  that  case  it  is  shown  in  Reference  26  that  only  about  8  percent  of  the 
resultant  jet  velocity  remains  effective  at  the  profile  position  equivalent 
to  the  boundary-layer  displacement  thickness.  Similar  factors  are  derived  in 
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Reference  17.  Thus,  on  a  somewhat  intuitive  basis,  the  incremental 
contribution  to  the  local  pressure  coefficient,  caused  by  the  effective 
"cavity"  velocities  over  the  bound  legs  of  the  vortex  feeding  sheet,  can 
be  computed  by  modifying  that  portion  from  Equation  (14)  as  follows: 


■2irpK  AL.  .C 
H  n.  .  in  w 


q  S.  . 


.  .Fv  . 
ci|  mj 


T 


:  V 
o  « 


(18) 


where:  q  is  the  dynamic  pressure;  K  is  the  final  augmented  circulation 

nij 

strength  of  the  n“  vortex  element  at  span  station  i,  and  AL. .  is  the 

corresponding  bound  vortex  link  length;  and  S.j  is  the  area  of  the  wing 

lattice  panel  affected  by  K  ...  The  factors  k  .  and  Gw  are  required  by  a 

V-L  representation  because  of  the  extreme  differentials  in  effective  outer/ 

inner  and  upper/lower-surface,  respectively,  flow  velocities.  The  k  . 

adjusts  for  the  effective  velocity  sensed  by  the  bound- lattice  network  in 

the  upper  surface  due  to  the  velocity  gradient  normal  to  the  surface  in  the 

cavity  backflow  region.  This  gradient  is  related  to  the  depth  of  the  cavity 

formation  region,  proportional  to  sin  a  ,  over  which  the  velocity  change 

occurs.  If  an  "effectiveness  factor"  is  defined,  based  on  Ref.  26  data,  as 

e  .08,  then 
0 


and 


(19) 


Cw  =  Average  velocity  over  the  upper  surface  ~  £ 


(20) 


in  the  case  of  augmented  lift  with  reattached  flow  over  much  of  the  upper 


surface. 
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The  factor  eQ  is  also  applied  in  determining  the  effectiveness  of  the 
wall  jet-like  profile  when  it  is  penetrating  the  forward  airfoil  region 


of  greater  backflow  velocity  due  to  the  vortex  tangential  component,  VQ 

^ni 

Tn_s  amounts  to  a  reduction  in  effective  free  stream  cross-flow  velocity, 
explaining  in  part  tne  concluding  results  in  Equations  (5)  and  (15)*  The 
available  empirical  formulas  only  allow  a  mean  cross-flow  velocity,  V  ,  to 
be  applied  to  the  jet  velocity  and  trajectory  computations.  However,  as 
snown  in  Appendix  A,  after  this  initial  approximate  pattern  is  established, 
experience  indicates  that  improved  values  of  Vg  can  be  obtained  at  each  span 
station  i  t?y  using  the  relationship 

ve  =  V®  cos  (0  +  A)  -  e  V  (21) 

i  ^ni 


across  the  wing  span.  The  vortex  tangential  velocity  at  a  distance 
(rni  +  Sj)  below  the  "locked"  or  suspended  leading-edge  vortex  axis  (refer 
to  Figures  6  and  9)  is  approximated  as, 


V  *  K  ./(r  .  +  8  )  +  Kt  ./r  . 
0  ■  nr  '  ni  Jy  Jnr  ni 
ni 


Kni  =  -2 ii 


kt  .  =  Ar_ ,  /0 

Jm  jhi/2  7 r 

as  the  jet  circulation  augmentation  from  equation  (12) 


Finally,  the  vortex  effective  radius  r^  is  defined  in  Appendix  A  along  with 

a  discussion  of  the  viscid  method  for  computing,  K  ,  the  final  circulation 

nij 

strength  variation. 
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3* 5 *3* 3  Jei-Vort£x_V^s£Oii8-S^abil_it7_Effe_c_ts  -  The  J-V  viscous  effects 
program  also  computes  vortex  burst  criteria  as  delineated  in  Appendices  A 
and  B.  For  wing  stations  past  the  critical  stability  points,  highly 
turbulent,  separated  flow  is  predicted.  Thus  the  values  of  C^,  lift  and  moment 
from  potential  theory  must  be  modified.  Fortunately,  scientists  at  the 
Loclcheed-Georgia  Research  Laboratory  devoloped  (in  Reference  12)  a  modified 
Helmholtz  theory  more  generally  applicable  to  stalled  lifting  surfaces. 

Basically,  the  method  employs  the  procedure  described  in  Reference  9»  but 
it  uses  a  modified  upper  surface  suction  pressure,  logically  deduced  from 
known  drag  characteristics  in  tne  presence  of  extensive  separated  flow  over 
airfoil  sections.  Similar  crossflow  drag  concepts  have  been  adopted  success¬ 
fully  for  nonlinear  lift  representation  in  References  32  and  33*  l'he  assumption 
is  made  that  the  revised  average  near-wake  pressure  applies  in  the  separated, 
turbulent  flow  region  created  by  vortex  burst.  Then  the  ratio  of  the 
Helmholtz  lift  to  the  potential  lift,  from  Reference  12,  is  as  follows, 

H  =  2.  sin(2o)/(4  +  it  sin  a)  sin  a  (25) 

8 

This  factor  is  multiplied  times  the  C  from  the  potential,  vortex-  lattice 

*V-L 

solution  as  influenced  by  the  appropriate  supervelocities.  At  high  angles  of 
attack  this  method  is  remarkably  successful  as  shown  in  Section  4*0» 

THEORETICAL  RESULTS  AND  COMPARISON  WITH  TESTS. 

At  lower  angles  of  attack,  when  vortex  burst  occurs  in  the  stall-transition 
region,  an  improved  correction  can  be  obtained  by  taking  the  ratio  of 
Helmholtz-to-potential-flow  lift  directly  from  Figure  4  at  the  given  angle. 

At  high  incidence,  when  the  jet  and  vortex  are  swept  aft  quite  severely 


because  of  too  little  mass-momentum  injection,  the  primary  vortex  may  be 
weak  and  too  small  to  prevent  a  large  separation  region  from  forming  (due  to 
secondary-vortex/turbulence)  forward  of  it  and  aft  of  the  leading-edge. 
Preliminary  calculations  at  =  .115  indicate  that  it  is  pertinent  to  apply 
the  modified  Helmholtz  factor  in  this  region,  as  well  as  outboard  of 
burst-critical  stations,  prior  to  the  addition  of  the  incremental  pressures, 

AC  ,  from  the  vortex  feeding  sheet  analysis. 
pv 

Por  some  cases  at  higher  incidence  and  C ^  ,  the  Owen's  time/spacially  varying 
eddy  viscosity  formulation  (Refs.  14  and  25),  combined  with  the  Rossby  number 
stability  criteria  (Appendix  A  and  Ref.  54)  seems  to  be  capable  of  predicting 
some  conservation  of  angular  momentum  after  burst  occurs.  As  the  aging 
process  re- initiates,  the  vortex  again  re-coalesceses  into  a  coherent 
rotating  fluid  mass.  If  this  occurs  on  the  wing,  it  is  possible  that  the  V-L 
computed  pressures  are  again  satisfactory  past  the  region  of  burst  recovery. 

Additional  confidence  is  gained  in  the  viscous  methods  employed  by  comparing 
tne  calculated  time-dependent  turbulent  eddy  viscosity  v  with  available 
test  data  and  with  recent,  more  sophisticated  methodology.  Figure  12  shows 
that  the  turbulent  eddy  viscosity  ratio,  vj/v  ,  predicted  herein  agrees 
favorably  with  wing  vortex  wake  data  summarized  in  Ref.  55  as  a  function  of 
vortex  Reynolds  Number,  \'/v.  The  analytical  point  is  for  the  8-degree  angle- 
of-attack  case  where  the  leading-edge  vortex  quickly  leaves  the  wing  upper 
surface,  and  penetrates  about  one  wing  span  into  the  near  wake,  just  before 
vortex  burst  occurs.  Similarly,  comparisons  were  made  with  wake  vortex  radius 
growth  predictions  using  boundary  layer  turbulent  shear-stress  concepts  by 


PRESENT  VISCOUS  VORTEX  METHOD 
(NEAR-WAKE  AT  MARGINAL  STABILITY) 


(a  ■  8°) 


TEST  DATA  RANGE 
(REF.  35  ) 


/ 


Dr.  J.  F.  Nash  (Conuultant  to  Lockheed,  Hof.  $6).  Again,  qualitatively 
good  agreement  wue  obtained  on  intermediate  wake  vortex  growth  rates  since 
both  methods  showed  similarly  significant  differences  from  the  classical 
"aquaro-root-of  :t"  variation,  depending  upon  model-scale  effects. 


4.0  THEORETICAL  RESULTS  AND  COMPARISON  WITH  TESTS 


4.1.0  LOW  MCHMCE  RESULTS  OF JTORTEX—LIET  CONTROL  THEORY 

A  distinct  difference  in  the  character  of  jet-vortex  interaction  occurs 
depending  upon  the  presence  of  leading-edge  separation.  The  mathematical 
formulation  was  primarily  based  upon  test  flow  visualization  (Refs.  2,  7 

and  8)  at  high  airfoil  incidence  where  separated  flow  exists  as  sketched  in 
Figure  5  (inset  view  for  stalled/Helmholtz  flow)  or  in  Figure  8.  However, 
at  low  incidence,  leading-edge  separation  is  either  non-existent  or  too 
insignificant  to  sustain  a  vortex.  As  discussed  in  Section  5*3  (ana  visualized 
in  Figures  30  and  31)  of  Volume  I  to  this  report,  even  at  low  angle  of  attack 
the  high  cross-flow  shear  layer  above  the  jet  propagates  vorticity  which 
causes  the  jot  region  to  develop  a  high  swirl  rate  as  it  is  swept  aft  over 
the  wing.  Note  that  the  analysis  is  limited  to  wing  root  blowing  cases. 

4*1.1  Cross-Flow  Jet  and  Vortex  Trajectories 

The  mean  trajectories  of  the  jet  and  vortex  (center  of  swirl)  were  estimated 
from  the  flow  visualization  photographs  cited  previously  and  plotted  in 
Figure  1J.  The  corresponding  jet  and  vortex  paths,  as  approximated  in 
Appendix  A,  are  shown  in  the  figure  to  be  in  very  good  agreement  with  the 
test  data  as  interpreted  for  C^  ~  .4  and  a  =  8.5  degrees.  Surprisingly,  the 
simplified  vorticity  integration  performed  herein  predicts  that  the  mean 
path  of  concentrated  vorticity  crosses  over  and  aft  of  the  jet  apparent 
trajectory  for  an  angle  of  attack  below  that  for  natural  airfoil  stall.  The 
analysis  makes  the  arbitrary  assumption  that  the  vortical  layer  extends  aft 
from  the  wing  apex,  or  the  junction  of  the  wing  leading-edge  and  fuselage 
sidewall.  This  is  not  apparent  in  the  test  results  except,  perhaps,  at  high 
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angles  of  attack  (c.f.  Volume  I,  Figure  27). 

4.1.2  Differential  Pressures 

At  low  incidence >  the  accuracy  of  the  basic  vortex-lattice  (V-L)  program 
is  examined  in  Figure  14(a).  At  =  0  and  a  =  0.5  degrees,  the  flow  is 
naturally  attached  and,  as  anticipated,  the  differential  pressures  from 
Volume  I  measurements  agree  quite  well  with  the  analysis.  The  leaning-edge 
singularity  is  cnaracteristic  of  all  lifting  surface  theories.  Although 
approximate  corrections  could  be  made  for  the  finite-radius  leading-edge  of  the 
wing  and  for  the  minor  thickness  or  camber  effect  observed  at  about  mid¬ 
chord  position  in  the  tests,  this  was  not  deemed  to  be  of  sufficient 
importance  for  the  objective  of  predicting  power-augmented  ?.ift  effects. 

For  a  spanwise  blowing  momentum  coefficient,  =  O.46,  the  comparison  of 
chordwise  differential  pressure  variations  at  36.4  percent  of  span,  as  in 
Figure  14(b),  indicates  that  the  theoretical  simulation  generally  predicts 
lift  increment  trends  like  those  measured  in  the  wind  tunnel  tests.  The 
effective  thickness/camber  increase  over  the  aft  chord  region,  under  the  jet/ 
vortex  path,  is  similar  in  distributed  form  to  that  from  the  tests,  but  is 
slightly  deficient  in  magnitude  and  locations  of  pressure  peaks.  Some 
of  this  may  be  due  to  the  relatively  coarse  spanwise  grid  spacing  (5  spanwise 
and  10  chordwise  collocation  stations)  used  in  the  lattice  (V-L)  program  to 
minimize  computer  time.  In  the  case  of  highly  curved  trajectories,  the  jet 
and  vortex  induced  velocities  were  diagonally  distributed  over  a  wing 
collocation  panel,  in  close  proximity  to  lattice  singular  points.  A 
smoothing  process  described  in  Section  3*5  minimized  the  possible  errors 
caused  by  such  a  "lumped-parameter"  system,  but  some  biased  weighting  is 


inevitable  in  the  correction  formulas  employed.  It  is  more  likely  that  the 
approximations  in  vortex  strength  and  jet  displacement  thickness  cause  the 
major  discrepancies  observed  in  this  non-optimum  case.  In  the  optimum  lift 
augmentation  cases,  with  a  highly  separated  flow  region  and  less  effective 
cross-flow  velocity,  the  error  in  calculated  jet  displacement  thickness  is 
probably  less  significant  to  the  overall  lift  predictions. 

4.2.0  HIGHjraCpMCEJffiSULTSJDP  VORTEX-LIPT  CONTROL  THEORY 

Leading-edge  or  primary  vortex-lift  augmentation  is  shown  in  Volume  I  to  be 
effective  in  increasing  maximum  lift  for  attitudes  above  normal  stall  and  for 
wing  root  blowing  momentum  coefficients  greater  than  about  0.1.  Lesser  amounts 
of  Cjj  do,  of  course,  provide  lift  augmentation  by  other  means,  as  discussed 
in  Section  4.1.0.  A  high  value  of  0. 46, above  the  optimum,  was  selected 
for  most  of  the  analytical  studies  since  it  was  desired  to  emphasize  the 
strong  vortex  formation  and  control  mechanisms. 

4.2.1  Cross-flow  Jet  and  Vortex  Trajectories 

The  experimental  jet  and  vortex  trajectories  are  presented  in  Figure  15,  as 

interpreted  from  oil  flow  photographs  in  Figure  20  of  Volume  I  at  a  slightly 

lower  value  of  C^  .  The  analytical  approximation  for  reducing  the  effective 

y 

cross-flow  velocities,  e^  across  the  span  yields  very  good  agreement  with 
the  test  observations  for  both  vortex  and  jet  paths.  The  jet  effective 
spread  angle,  y  =  20  degrees  used  in  the  analysis  is  not  directly  determinable 
from  the  surface  flow  pictures.  But,  as  pointed  out  in  Volume  I,  if  the  jet 
direction  and  spread  rate  is  not  proper,  the  vortex  will  not  be  positioned 
properly  behind  the  wing  leading-edge.  For  the  oval  nozzle  used  in  these 
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tests  (see  Volume  I),  the  geometric  exit  angle  of  12  degrees  aft  relative 
to  the  wing  leading-edge  sweep  angle  was  not  the  observed  initial  jet 
trajectory  angle,  #Q,in  the  photographs.  Further  study  is  needed  on  the 
effects  of  nozzle  wall  and  orifice  unsymmetries.  For  this  analysis,  the 
initial  angle  was  adjusted  to  match  that  in  the  photographs  in  order  to 
ensure  reasonable  starting  conditions.  Fortunately,  no  other  adjustments  in 
coefficients  for  the  semi-empirical  trajectory  solutions  presented  in 
Appendix  A  were  required  to  achieve  the  good  agreement  in  Figures  15  and  15 
for  vortex  location  versus  span  over  the  range  8,5°<.  a  <  27°  and  0.1<  C^<  O.46. 

4. 2. 2.0  Stall/Burst/Reynolds  Number  Effects 

Vortex  burst  was  predicted,  accompanied  by  a  rapid  increase  in  apparent 
vortex  radius,  at  about  80  percent  of  span,  using  the  Rossby  number 
criteria  (Appendix  A).  The  oil  streak  lines  in  Figure  50  of  Volume  I  also  show 
a  sudden  divergence  in  this  outer  span  region.  Simultaneously,  the  pitching 
moment  for  similar  and  angle  of  attack  (27°)  in  Figure  29  of  Volume  I 
shows  a  sudden  reversal  in  trend,  toward  nose-up  conditions.  Thus  the  tests 
indicate  the  possibility  of  partial  wing-tip  stall  initiated  by  vortex  burst. 

This  same  case  was  later  rerun  with  a  factor  of  twelve  increase  in  wing-scale, 
and  the  Rossby  number  did  not  predict  vortex  burst.  Yet  the  differential 
pressures,  prior  to  modification  for  burst/separated  flow  at  small-scale,  as 
described  in  Section  5*5«3>  showed  only  minor  changes  with  scale.  Professor 
Owen's  dissipation  analysis  employed  herein  (Appendix  A)  does  predict  changes 
with  scale  in  the  vortex  ageing  process.  This  predicted  improvement  in 
spanwise  blowing  at  large  scale  is  partially  substantiated  by  larger  scale 
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model  tests  discussed  in  Section  5*2.1 »  Volume  T,  since  slightly  less 
momentum  coefficient  seems  to  be  required  for  a  given  lift  increment. 
Unfortunately,  some  differences  in  configuration  existed  in  these  Lockheed- 
funded  tests,  so  no  firm  conclusions  can  be  drawn. 

Some  discussion  of  the  effects  of  the  natural  stall  angle  for  a  given  air¬ 
foil  (in  this  case,  14  degrees),  on  the  two-  and  three-dimensional  analyses 
has  already  been  given  in  Sections  2.2  and  3.5*  Also,  the  starting-vortex 
strength  limitations  of  Figures  5-7  must  be  kept  in  mind.  Then  the  vortex- 
lift,  ^CLV,  variations  with  angle  of  attack  for  two  extremes  of  (0.115  and 
O.46)  can  be  calculated  from  the  3-D  simulations,  pre  lented  as  dashed  lines  in 
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Figure  16.  The  "unstalled"  limit  angle  is  again  shown;  but,  in  addition,  the 
vortex-burst  initiated,  partial  wing  stall  region  is  approximately  indicated. 
This  effect  occurs  at  high  incidence  where  the  is  only  marginally  sufficient 
to  entrain  fluid  and  propagate  fully  through  the  leading-edge  separation 
region,  causing  local  chordwise  reattachment.  The  test  points  in  Figure  16 
at  similar  1  e  (0.155  and  O.46)  show  the  tendency,  as  hypothesized  in  the 


analysis,  to  reach  an  asymptotic  ACT  limit  as  this  burst/stall  condition 

L 
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is  approached.  At  lower  angles  of  attack,  the  differences  relative  to  theory 
are  a  result  of  the  difficulties  in  separating  out,  from  the  test  data,  the 
lift  increment  due  to  the  vortex  alone,  as  estimated  for  Figure  32  in  Volume  I 
of  this  report.  This  problem  is  even  more  complex  when  combined  inboard  and 
outboard  spanwise-blowing  jets  are  used. 


4* 2. 2.1  Lift  Increments  and  Causes  -  The  total  lift  from  potential,  vortex 
lattice  theory  (Refs.  22-24)  compares  well  with  the  test  data  prior  to  stall, 
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FIGURE  16  COMPARISON  OF  BOTH  LOW-  AND  H I GH- 1 NC I DENCE  STALL  FOR  ACL  IN 
TESTS  WITH  CORRESPONDING  EFFECTS  FROM  THE  ANALYSIS  v 
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as  shown  in  Figure  17(a)  at  «=  0.  This  does  imply  compensating  effects, 
concerning  the  natural  wing  camber  and  thickness  corrections  for  this 
Lockheed  super-critical  airfoil,  for  middle  range  fit's,  since  the  V-L  method 
does  not  include  these.  A  non-planar  wake  is  considered,  however,  giving 
some  nonlinear  induction  in  the  lift. 

Obviously,  any  powered-lift  technique  which  reattaches  the  flow  at  above-stall 
angles  will  yield  significant  AC.  increments,  shown  as  "Indirect  Augmentation" 
in  Figure  17(b).  This  probably  consists  largely  of  a  boundary  layer  control 
effect  related  to  the  spanwise  blowing  jet  spreading  over  the  upper  surface 
aft  of  the  reattachment  zone.  Essentially,  this  re-establishes  the  potential 
flow  character,  with  the  addition  of  some  effective  camber  due  to  the  jet 
displacement  thickness, shown  by  the  dotted  line  in  Figure  17(b)*  Refer  to 
the  upward,  nearly  parallel  shifts  in  the  versus  a  curves,  at  zero  incidence, 
in  Figure  29  of  Volume  I. 

The  increase  in  the  slope  of  the  lift  curve  above  the  stall  is  directly 
attributed  to  the  vortex  and  its  induction  effects  through  the  elements  of 
the  feeding  sheet  (dashed  line  in  Figure  17(b)).  Actually,  most  of  the 
lift  above  the  mean  stall  level  exists  due  to  the  reattachment  caused  by 
the  formation  of  at  least  a  small  vortex  near  the  leading-edge,  so  the  above 
iB  an  arbitrary  division  of  augmentation  increments  due  to  the  vortex. 

At  very  high  angles  of  attack,  the  vortex-burst-created  stall  over  part  of 
the  span  causes  some  lift  loss  (dot-dashed  line,  Figure  17(b)),  and  this 
agrees  quite  well  with  the  test  stall  character  above  a=  25  degrees  (solid 
line  in  the  figure). 
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FIGURE  17 


COMPARISON  OF  INCREMENTS  OF  LIFT  FROM  VORTEX-LIFT  AUGMENTATION  ANI 
SPANWISE  BLOWING  IN  BOTH  TESTS  AND  ANALYSES 


4. 2. 2. 2.  Differential  Pressures  -  The  previously  discussed  wing  tip  stall 
phenomenon  is  revealed  m  part  by  the  non-potential  flow  character  of  the 
aft-chord  distribution,  at  'Jl  percent  of  semi-span,  shown  in  the  top  part 
of  Figure  IB.  The  nearly  flat  aft  distribution  is  more  like  that  for  Helmnoitz 
flow,  used  to  modify  the  V-L  simulation  for  wing  regions  encompassing  the 
indication  of  vortex  burst  (See  Section  3.5.5). 

At  72.7  percent cf  semi-span,  the  aft-chord  differential  pressures  are  again 
i.n  agreement  between  test  and  theory.  There  is  some  indication  that  the 
simulation  under-predicts  the  leading-edge  vortex  strength  at  this  station; 
but,  unfortunately,  at  critical  chordwise  stations  the  pressure  measurements 
were  consistently  unreliable,  as  indicated  in  the  figure. 

At  36.4  percent  of  semi-span,  the  vortex  induced  peak  1b  predicted  quite 
accurately;  and  tnese  C^'s  probably  would  improve  if  the  wing  leading-edge 
singularity,  from  the  V-L  method,  were  corrected  for  a  finite  leading-edge 
radius.  Other  minor  differences  may  also  improve  if  a  greater  number  of 
collocation  panels  are  selected,  as  discussed  relative  to  Figure  14  (Section  4»1)> 

Similar  agreement  between  test  and  analysis  is  shown  in  Figure  19  for  an 
intermediate  angle  of  attack,  a=  21  degrees,  and  =  0,44*  At  this  outboard 
span  station  the  simulation  seems  to  over-predict  the  jet-cross-flow  super 
velocity  effect  over  aft  chord  sections,  and  a  smoother  transition  between 
the  vortex  and  jet-swirl  pressure  peaks  is  indicated  than  that  measured  in  the 
wind  tunnel  model  tests. 

The  pressure  distributions  in  Figures  14  and  18  reveal  some  interesting  fea¬ 
tures  of  the  flow  over  the  upper  surface  of  the  wing.  For  spanwise  blowing 
at  low  wing  incidence,  the  aft-chord  jet  swirl  peak  is  very  apparent.  But, 
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FIGURE  18  COMPARISON  OF  CHORDWISE  PRESSURE  DISTRIBUTION  FOR  THEORY 
AND  EXPERIMENT  WITH  SPANWISE  BLOWING  OVER  A  SWEPT  WING 
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FIGURE  19  COMPARISON  OF  CHORDWISE  PRESSURE  DISTRIBUTION  FOR  THEORY  AND 
EXPERIMENT  WITH  SPANWISE  BLOWING  OVER  A  SWEPT  WING 
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at  high  attitudes,  the  influence  on  the  pressure  coefficient  caused  by  the 
vortex  and  feeding  sheet  near  the  leading-edge  is  much  more  dominant.  The 
increased  suction  is  partly  attributable  to  the  higher  upwash  velocities 
forward  of  and  at  the  vortex  center  within  the  "cavity"  or  decelerated  mean 
flow  region.  As  the  reattachment  point  is  approached,  aft  of  the  vortex  for 
tne  most  inboard  station  shown  in  Figure  18,  the  pressure  increases  rapidly 
due  to  downwash  velocity  components  from  both  the  vortex  and  jet  impingement 
on  tne  surface  forward  of  the  nozzle.  However,  the  jet  wake  provides  a 
cross-flowing  super-velocity  to  elements  of  the  upper  surface  vorticity 
feeding  sheet  which  then  causes  anotner,  smaller,  suction  peak  to  form 
xartner  aft  on  the  chord.  These  effects  are  less  pronounced,  but  are  still 
apparent  at  more  outboard  stations.  The  finite  displacement  thickness  of 
tne  jet  suggests  a  local  camber  effect;  but  the  AC^'s  are  higher  than  could 
be  expected  purely  by  this  meclianism. 


The  three-dimensional  simulation  uses  a  feeding  sheet  with  vortex  elements 
wruch  can  allow  the  jet  to  produce  pressure  variations  qualitatively  consistent 
with  the  above  experimental  observations.  Since  quantitative  agreement  of  the 
modeling  technique  witr.  Polhamus'  leading-edge  suction  analogy  for  delta  wings 
(Refs.  2o  and  21)  has  been  shown  (Figure  10),  future  research  would  be 
desirable  to  determine  if  some  combination  of  these  methodologies  could  be 
used  to  predict  the  pressure  distributions  (and  pitching  moments)  qua  litatively 
for  more  general  planforms  either  v/itn  or  without  powered  lift  rfiects. 


4.3.0  SPANWISE  LIFT  DISTRIBUTIONS  AND  TOTAL  LIFT 


Tne  changes  in  spanwise  lift  distribution  caused  by  vortex-lift  augmentation 
are  apparent  m  Figure  20.  The  potential  (V-L)  theory  at  high  incidence, 
a=  2 /°,  arbitrarily  assumes  attached  flow  to  give  a  classical,  near- 


elliptic  distribution.  As  normalized  to  the  Maximum  in  each  case,  and 
using  a  blowing  momentum  coefficient,  C ^  =  O.46,  the  outboard  lift  drops 
relative  to  that  at  inboard  sections.  Part  of  this  is  due  to  wing-tip 

Q 

stall,  as  previously  discussed.  But  the  calculated  inboard  C*—- 

c 

peak  is  partly  due  to  the  vortex  augmentation  and  partly  due  to  thrust 

recovery.  For  example,  an  integration  of  the  test  pressure  data  at  36.4 

0 

percent  of  semi-span  yields  only  the  normal  force  coefficient  tr,  and 
the  inboard  hump  is  not  apparent  in  this  value.  The  leading-edge  thrust  and 
the  drag  at  that  wing  station  are  difficult  to  obtain  because  of  the  few 
pressure  taps  available.  Thus  a  rough  calculation  was  made  for  the  lift- 
vector  components  of  thrust  and  drag  from  the  theory,  and  then  combined  with 
the  vertical-component  of  test  normal  force  at  a  =  27°.  This  estimated 
test  point  is  shown  in  Figure  20  as  a  dotted  Cxrcle.  Relatively  good  agreement 
oe tween  test  and  analysis  is  then  indicated.  However,  the  fuselage  correction 
in  the  simulation  is  very  approximate  (See  Section  3*!?*3)»  so  the  most  inboard 
part  of  the  carve  should  not  be  used  quantitatively.  Intuitively,  a  decline 

Q 

of  —  in  this  region  would  be  expected  at  these  high  angles  of  attack, 
c 

Comparison  oi  test  and  analysis  for  C.  versuB  C „  in  Figure  21  shows  very 

LMAX 

good  agreement.  Empirical  adjustments  to  the  two-dimensional  starting 
conditions  should  be  recalled,  however,  so  as  to  not  give  misleading  con¬ 
clusions  on  the  state-of-the-art  in  this  vortex-lift  augmentation  theory.  Also, 
it  is  fortunate  that  the  classical  camber  and  thickness  corrections  apply  so 
well  to  these  test  cunditions.  For  other  model  configurations,  the  modified 
vortex-lattice  program  (Figure  11)  has  no  basic  limitations.  But  a  more 
analytical  formulation  of  the  cross-flowing  jet  and  v: scous/boundary-layer 
affects  is  still  necessary  for  more  general  applications. 
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FIGURE  20  COMPARISON  OF  LIFT  DISTRIBUTION  FROM  PRESENT  THEORY 
WITH  THAT  FROM  POTENTIAL  THEORY  AT  =  0 


TEST  (OVAL  NOZZLE.*  INBOARD 
\  MOMENTUM  RATE 
\  INJECTION) 


dQ  =  .28  IN. 


/r  V THEORY  (VORTEX- LATTICE  WITH  VORTEX  BURST 
"  HYPOTHESIS  FOR  PARTIAL  STALL) 


NOTE:  DATA  INCLUDES  CLASSICAL  WING  AND  JET 
EFFECTIVE  THICKNESS  AND  CAMBER 
CORRECTIONS  TO  ALLOW  COMPARISONS. 
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FIGURE  21  EFFECT  OF  SPANWISE  BLOWING  JET  ON  Clmax 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 

5 . 1  CONCLUUINGJffiMABJ^ 

The  three-dimensional  analysis  for  vortex-lift  control,  presented  herein, 
is  admittedly  very  empirical  in  nature.  However,  the  presence  of  each  of 
the  hypothesized  elements  of  the  phenomena  inherent  in  spanwise  blowing  is 
strongly  indicated  by  the  experimental  evidence  delineated  in  Volume  I  of 
this  report, 

A  previous  analytical  solution  for  viscous  vortex  flows,  verified  by  idealized 
experiments,  also  has  been  presented  in  revised  form  so  as  to  show  character¬ 
istic  dependence  on  radial-inflow  Reynolds  Number  and  Rossby-number  concepts. 
This  describes  a  fundamental  physical  mechanism  for  causing  an  organized  vortex 
to  grow  from  an  upstream,  separated  shear-layer  and  to  be  stabilized.  In 
addition,  if  a  separated  flow  region  exists  or  can  be  created  adjacent  to  this 
feeding  sheet,  then  adverse  pressure  gradients  may  prevent  any  major  (periodic 
or  steady)  streamwise  shedding  from  a  predominantly-spanwise  vortex.  Such 
adverse  pressures  may  be  created,  for  instance,  by  aspanwise  jet  just  down¬ 
stream  of  the  vortex  formation  region.  Over  a  wing  upper  surface,  this 
results  in  a  strongly  sheared  three-dimensional  wall-jet  velocity  profile. 

An  appropriate  Stokes'  integration  of  this  wall-jet  vorticity,  in  planes 
oriented  so  as  to  yield  increased  circulation  to  the  primary  vortex,  provides 
the  correct  amount  of  vortex-lift  augmentation  in  comparison  with  experimental 
data.  The  ensuing  viscous  dissipation,  calculated  using  a  realistic  and  time- 
dependent  turbulent  eddy  viscosity,  does  not  significantly  reduce  this  lift 
augmentation  at  high  angles  of  attack. 
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The  prediction  method  for  vortex-lift  augmentation  derived  herein  requires 
initially  a  two-dimensional  source-sini/vortex  model  for  locating  the  leading- 
edge  vortex.  In  addition,  a  first  estimate  is  made  for  this  vortex  circulation 
strength  consistent  with  the  appropriate  boundary  conditions  and  momentum 
coefficients.  Then,  with  available,  approximate  cros3-flowing-jet  solutions 
for  trajectory  and  entrainment,  the  spanwise  distribution  of  the  equivalent 
two-dimensional  singularities  are  estimated.  These  induction  effects  are 
combined  with  those  from  a  spatially  varying,  finite-element  vortex  feeding 
sheet  emanating  from  the  wing  leading  edge.  Subsequently,  the  combined 
supervelocities  are  employed  in  a  vortex-lattice  solution  for  general  wing 


configurations.  The  total  forces  and  aerodynamic  coefficients  are  compared 


witn  the  corresponding  test  data  from  Volume  I,  and  relatively  good  agreement 


is  obtained.  Optimum  mass-momentum  injection  conditions  can  be  predicted 


using  this  method,  for  initial  design  estimates.  Pa-vorable  Reynolds  Number 


effects  are  also  possible  since  a  large  scale  model  showed  less  tendency 


for  the  augmented  vortex  to  burst,  and  corresponding  flow  separation  over 


sections  of  the  wing  is  less  likely  to  occur.  Otherwise,  lift  coefficients 


predicted  for  small  scale  models  are  substantially  the  same  as  those  for  large 


scale  models,  similar  to  the  conclusions  from  scaled  tests  in  Volume  I  of  this 


report. 


The  combination  of  nonlinear  viscid  and  inviscid  mathematical  elements  used 


in  this  initial  formulation  require  that  an  iterative  and  semi-empirical 


procedure  be  employed  for  solution.  Now  that  the  heuristic  objectives  of  this 


research  program  have  been  satisfied  (i.e.,  to  obtain  an  improved  physical 


understanding  01  the  numerous  flow  mechanisms  in  spanwise-blowing/vortex  lift 


augmentation),  more  sophisticated  viscous  flow  theories  are  required  to 


reduce  the  degree  of  empiricism  and  thereby  to  enable  application  of  the 

above  methods  more  widely  and  with  {jreater  confidence. 

5.2  RECOMfffiOATIONS. 

The  following  tasks  are  recommended  for  future  research  to  advance  current 

understanding  of  the  spanwise-blowir.g  phenomenon: 

1.  Develop  a  rational  analysis  of  jet  flow  in  a  cross-stream  witn  a 
parallel  wall.  This  should  involve  both  mathematical  and  experimental 
investigation  of  jet  structure,  to  incorporate  measured  turbulence 
parameters  in  an  analytical  framework  for  predicting  jet  spread  rate 
and  decay  using  three-dimensional  boundary  layer  methodology. 

2.  Further  refinement  of  the  two-dimensional  theory,  incorporating  a  leading- 
edge  feeding  sheet  for  the  vortex  and  inclusion  of  rigorous  criteria 

for  selecting  vortex  position  uniquely  for  each  flow  condition. 

3.  Basic  research  is  needed  on  the  phenomenon  of  vortex  bursting, 
especially  in  the  presence  of  a  lifting  surface. 

4.  Finally,  the  results  of  items  (l)  -  (3)  should  be  used  to  refine  the 
present  vortex-lattice  scheme  for  three-dimensional  loads  to  the 
point  that  optimum  design  studies  can  be  made. 
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APPENDICES 


APPENDIX  A:  Summary  of  Approximate  JetAortex  Viscous  Flow  Methods. 

The  physical  problem  discussed  in  this  report  poses  so  many  complex  flow 
interactions  that  theoretical  analyses  of  even  sub-elements  of  the  formulation 
(as  outlined  in  Figure  11  of  Section  3.0)  have  only  been  accomplished  in  part, 
or  semi-empiricaily.  The  purpose  of  this  Appendix  is  to  summarize  approximations 
obtained  from  available  literature  and  employed  witn  minor  modifications  elsewnere 
in  tne  text. 

General  Fojmulas_£or_j£ts.  and  Vortices 

The  free-jet  in  a  cross  flow  is  still  usually  treated  semi-empirically,  although 
numerical  or  finite-element  methods  (e.g. ,  Refs.  15  and  18)  are  recently  mucn 
improved.  Unfortunately,  these  methods  are  cumbersome  and  lacking  in  consider¬ 
ation  of  the  boundary- layer,  viscous  wall-jet  interaction  most  essential  to  this 
problem.  On  this  basis,  the  momentum  entrainment  theory  developed  by  McAllister 
(Ref.  16)  seems  most  appropriate,  yet  simple  in  form.  Using  the  momentum 
theorem  and  specific  mass-momentum  flux  data,  he  showed  that  the  entrained 
mass  flux,  <q>,  carries  tne  undisturbed  cross-flow  velocity,  into  the  jet  in 
such  a  way  that  it  adds  the  momentum,  Q^,  to  the  jet  flow.  A  principal 
modification  employed  herein  is  the  replacement  of  by  V  ,  the  effective 
cross-flow  velocity  (Equation  21  in  Section  3.5),  defined  as 

Ve  *  Vos  (eo  +  A)  "  e0  \ 
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where 

#o  is  the  initial  jet  trajectory  angle,  as  shown  in  Figure  A. 1, 

v  is  the  initial  tangential  velocity  of  the  first  unbound  (leading- 
edge)  vortex  element,  at  a  height,  (y  -  r^),  above  tae  wing 
surface, 

A  is  the  wing  sweep  angle,  ana  eQ  is  a  surface-induction  effective¬ 
ness  parameter,  as  discussed  in  Section  in  this  Volume. 

(Refer  to  the  Nomenclature  Section  for  quantities  in  the  text 
but  not  redefined  here.) 

For  starting  conditions  in  the  wing  root  region  it  is  necessary  to  estimate 

the  value  of  equivalent  vortex  radius,  r^ ,  corresponding  to  the  initial  two- 

dimension  circulation  strength,  K  .  The  data  in  Figure  11  of  Volume  I  shows 

e 

that,  past  the  underexpanded  jet  length,  x  (determined  as  Y  from  Figure  13, 

s  n 

Volume  i)  the  spanwise  jet  decays  much  the  same  as  a  fully  expanded  jet,  for 

the  same  pressure  ratio,  ?S/p  •  If  the  initial  axial  coordinate  for  the  jet, 

am 

is  defined  in  terms  of  wing  coordinates  as  in  Figure  A.l,  then, 

h  =  xi/xm  cos  (A+  6o )  “  Vx  U*2) 

'  m 

where  x^  is  the  wing  semi-span  length  (o;b/2).  For  consistence  with  the 
terminology  of  Reference  16,  the  jet  "set-back"  factor,  a  s  0.8,  for  a  virtual 
origin  is  assumed  to  apply.  Also,  the  cross-flow  velocity  ratio,  redefined  as 
R  =  Vj/Vg, is  used  to  obtain  a  best-fitted  relationship  for  the  jet  spreading 
factor  Kj(r)  from  test  data  in  Reference  16  as 

Kj(R)  =  J  +  4//r 


(A.3) 


ThenMcAllister's  equation  for  jet  wake  depth  can  he  written  as, 


h.  as  (2aD  +  A  £  x  )Kt(R) 
i  '  o  ox  my  / 


(A.4) 


This  is  modified  in  two  ways.  Dq  is  the  effective  jet  diameter  at  the  beginning 
of  the  fully  expanded  region  given  by 

1_ 

I)  ^(P„/P  )2d  /. 

o  v  E'  am/  o  (A. 5) 

where  d  is  the  true  orifice  diameter  with  P_  and  P  as  the  jet  exit  and 
o  E  am  ° 

ambient  fluid  densities,  respectively.  Also,  Aq  is  an  empirical  growth 
coefficient  for  a  wall  jet  which  reattaches  over  a  curved  surface.  Prom  the 
LDV  measurements  (Volume  i),  and  some  data  in  Reference  17,  Aq  cs  0.04  is  a 
reasonable  estimate  for  present  purposes.  Since  the  jet  depth  tends  to  size 
the  leading-edge  "cavity"  region,  a  first  approximation  is 


n  a*  h±/2  (A.  6) 

for  £.>1, 

1 

so  that,  r1  a*  aDQ  (A. 7) 

This  is  not  critical,  since  the  n  are  recalculated  more  accurately  later. 


A  good  curve-fit  of  data  on  jet  trajectories  in  Reference  16  can  be  obtained 
by  the  expressions, 


zi  =  zj  +  xssin(0o+A)  +  7 rx 
+  xi  tan(@o+A) 


m 


1  -  cos1//2(^  (±) 


and 


(A.  8) 


6^  =  Arctan(dzi/dxi) 


(A.9) 
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where 

z.  is  the  streamwise  deflection  of  the  jet  relative  to  the  nozzle 

i 

.  exit  at  distance  zT  from  the  wing  origin,  at  wing  station  i 
and 

Q is  the  jet  turning  angle  at  the  same  station  i. 

Similar  adjustment  of  the  expression  (Ref.  16)  for  the  local  maximum  jet 

velocity,  B  ,  results  in 
i 


mi  6  D*s 

UJ  -  Ve  «i  + 

(A. 10) 

where 

D  =0'(1) 

(A. 11 ) 

for  the  cases  studied. 

The  axially  propagating  velocity  profile,  as  derived  by  Gortler  (Ref.  57) 
for  a  free  jet,  is  assumed  to  be  appropriate,  since  it  is  based  on  the  hypothesis 
that  the  eddy  viscosity  is  constant  along  a  crose-section  of  the  flow,  and 
this  is  consistent  with  initial  assumptions  on  viscous  effects  in  vortex 
flows.  The  equation  is  of  the  form; 


for  ^>0 

where  tne  argument  for  the  sech  is  modified  to  include  the  additional  jet 


spreading,  while  turning  over  a  surface,  observed  in  flow  visualization  tests 
(Volume  i).  Ag  —  r/ ^  appears  to  be  a  good  approximation  for  these  cases. 
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Additional  equations  are  adopted  (Ref.  16)  with  minor  modifications.  The 


volumetric  mass  flux*  Q^,  is 


^i  ~  1  +  |  J  cos  6  if  tan#.  -  tan  Q  )  +  ^iXi 

%  v  e  °/\  1  7  B  D, 


where 


and 


B  ~  9  “ 


&/x/9  + 


1) 


*  =  d  UT 
o  -  o  J 
4 


The  momentum  flux  is  approximated  as 


(A.13) 


(A. 14) 


(A. 15) 


where 


^i  _  cos  dp  -  io  Cf  lnpi  +  ^ 
Jq  ~  cos  $±  1 


J  ~  P  U.<i 
o  J  Jo 


(A. 16) 


(A.17) 

and  the  last  term  in  equation  (A.16)  is  based  on  a  similar  expression  by 
Hewman  (Ref.  17)  to  account  for  a  finite  friction  coefficient,  C^,  for  a 
wall  jet.  Equation  (A.17),  in  particular,  can  be  improved  upon  by  using  the 
gas  laws  for  compressible  flow;  but  these  expressions  were  found  to  be  useful 
for  first  estimates. 

For  vortex  dissipation  effects,  an  initial  estimate  of  axial  flow  velocity, 


in  outer-core  regions,  is  approximated  as  follows: 


vi 


Vx  .  -  V .  cos  e0  -  Ve  sin  °o 
vi  1l  1 


(A. 18) 


Hote  tnat,  from  this,  a  leading-edge  vortex  trajectory  angle  Q  can  be 


estimated, 
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/zT  z.  A  -1  [V  .  cos  (/, 

0V1  —  A  +  Qq  -  y /  +  tan  [ex 


'*V/vx  ‘ 

viJ 


(A. 19) 


Usxng  Equation  (22)  of  Section  3*5  in  Volume  II,  a  first  order  correction  to 


tne  vortex  augmented  strength,  Kvi,  is  calculated  as, 


K  .  ^ 
vi 


V  !  *ve/  \  Kv 

eu /  /ri  +  sj>J  6 


(A. 20) 


Then  Owen's  (Ref.  14)  aging  parameter,  Ai(t),  can  be  approximated  (in  the 
range  of  interest)  by  the  expression, 


where 


*  £  log1Q  (10  t ±/t) 


k=  i 

X. 


T  cos  "(XT#'"] 

vk  1 


(A.21) 


(A. 22) 


r  -  Xr,  (r,  +  8j)/Kv1 


(A.  23) 


The  widtn  of  the  turbulent  mixing  or  entrainment  region  8^  as  sketched  in 

Figure  (A. 2)  presents  a  complex  problem.  The  turbulence  in  this  rate-of-strain 

environment  must  increase  as  the  product  of  circumferential  velocity  VQ  and  r. 

e±  1 

increases,  according  to  Rayleigh's  stability  criterion  (Ref.  38) »  thus  causing 

tne  radius  and  the  mixing  region  to  grow  outward  into  the  more  stable  region 

where  the  product  •  r.  tends  to  decrease.  To  avoid  an  "eddy  viscosity" 

6i  1 

which  is  variable  with  scale  (radius  of  curvature),  Owen  (Ref.  14)  defines 


a .  =  6  •  /  r. 
1  1  /  x 


(A. 24) 


FIGURE  A-2  VORTEX  DECAY  MODEL  (REF.  14) 


and 


A  i  ii/ 


(A. 25) 


where  /  -s  the  local  average,  mixing  length  which  varies  along  the  vortex  axis. 
Trie  primary  vortex  strength  (Equation  A. 20)  is  augmented  by  the  spreading  jet 
front  as  it  advances,  the  net  effect  of  which  can  be  simulated  as  described 
in  Section  30  w^tn  tne  vortex  lattice  program.  However,  simultaneously  the 
turbulence  develops  a  Reynolds  shearing  stress  which  modifies  the  mean 
(irrotational)  motion,  generating  vortxcity  and  allowing  the  extraction  of 
kinetic  energy.  If  the  relation  (Ref.  14)  between  the  rate  of  dissipation  and 
tne  turbulent  velocity  q,  is  accepted  as: 


*SL 

dt 


-.1 

a 

tz 

i 

1 

(A. 26) 


where  k  is  a  proportionality  constant  of  order  unity,  integration  of  the 
equation  for  the  rate  of  change  of  turbulent  energy  With  respect  to  r  results  in 
an  expression  for  tne  turbulent  energy  balance  at  section  i  as  given  in 
Reference  14.  Subsequently,  it  is  shown  that  for 


T.  1 
TT  =  *i2  ’ 


(A.27) 


The  corrected  inner  core  radius  r  ^  (=r^in  Figure  A. 2)  is  given  as 
rc.  -  2\SV*i 


(A. 28) 


and 


6.  =  cr.  r 
1  xc. 

1 


(A. 29) 


Then  the  circulation  strength,  corrected  for  vortex  aging,  is  approximately, 
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K  .  : 
m 


,\3//nV\2  Ai(ti) 


(A.  30) 


For  an  efficient  spanwise  blowing  case,  the  vortex  radial  growth,  is  quickly 
stabilized.  For  turbulence  possessing  an  unobtrusive  radial  component  of 
velocity,  Professor  Owen  hypothesizes  that  k  and  X  «  5*  Since 


k  -  k  k 

i 

then,  the  scale  of  turbulence  is  about,  |.  s  5  ^  ,  indicating  that  the 

X  0 

1 

turbulent  eddies  may  be  predominantly  axial  in  character  in  the  mixing  region. 
This  seems  reasonable  in  view  of  the  high  axial  shear  flow  created  in  outer 
vortex  layers  by  the  coaxial  jet  in  the  case  of  spanwise  blowing. 


The  point  souice-sink  representation  for  velocity  components  is  of  standard  form, 
taken  directly  from  Reference  39 »  so  further  elaboration  here  is  unnecessary. 
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numerous  vortex  burst  nypotneses  or  criteria  have  been  studied,  since  it  was 
anticipated  teat  suen  instabilities  would  have  a  significant  effect  on  the 
^red.cted  vortex  lift  augmentation.  The  most  consistent  criterion  was  found 
to  be  tne  Rossby  number,  or  the  ratio  of  linear  ( inflow )~to-angular  momentum 
in  tne  vortex  (References  1,  34  and  4L>).  Reference  34  gives  a  particularly 
good  discuss xon  of  this  parameter  (more  generally  used  in  meteorological 
studies);  ana,  for  a  great  variety  of  flow  conditions  the  critical  level  _s 
snown  tc  be, 


hQ  =  0.261  (a. 32) 

wnere  a  distinct  change  mi  tne  viscid/mviscid  character  oi  the  flow  occurs. 

1ms  .i  magnitude  is  quite  well  verified  by  experiments  ir.  References  34  and  40. 

altaougn  for  complex  flow  situations  the  matnematical  expression  is  subject  to 
interpretation,  for  present  purposes  tne  iioesoy  number  is  calculated  as  a 
function  of  position,  i,  along  the  vortex  axis  as  follows: 

%  =  jVVfe  (A.  33) 

1  2jrrc  X,Kni 
1 

This  was  found  to  be  quite  successful  in  predicting  vortex  burst  as  observed  in 
the  test  data,  as  discussed  in  Sections  3.3  and  4.2  of  this  Volume. 


Another  criterion,  by  Lamboume  and  Bryer  (Ref.  41),  tended  to  predict  burst 
at  a  similar  wing  station  using  tne  formula 


V.  -  V 
Xv.  Xv(i+1) 

t. 


absolute 


(A. 34) 


1 


i 

% 
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Unfortunately,  multiple  bursts  were  sometimes  predicted,  at  unlixely  Wj.ng 
stations,  Ine  problem  may  arise  from  the  difficulty  _n  selecting  the  proper 
"stream- tube"  radius,  by  this  concept,  for  calculating  V^.  Also,  many  criteria 
are  based  on  criticality  concepts  applicable  to  other  tnan  burst  conditions, 
as  discussed  aptly  by  Kali  (Ref.  42). 

3ume  scale-effects  were  predicted  using  the  Radial  (inflow)  Reynolds  N  tracer,  as 


Q.-Q 

N  =  1  0 

1  i  27IVX. 

1 


(A. 35) 


*i;  is  defined  as  negative  for  inflow. 


finally,  the  vortex  Reynolds  Number  is  given  by 


K  'ty*  • 


U- 3<>) 
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iiPPi2.1I/.  B:  Background  Analysis  for  Wake  Vortex  Stabil.zation  Concept 


Reference  1  presents  a  tneoretical,  viscous  analyses  of  the  havier-otoKes 
equations  ior  an  axisymmetrJ.c  vortex  forming  in  a  near-wake  or  separated  flow 
refaJ.on.  ±>y  means  oi  a  transformation  oi'  tne  equations  and  a  vortex-growtn 
similarity  concept,  nonlinear  solutions  are  found  m  terms  of  vortex 
staDi.ii.ty  modes.  In  particular,  tne  Strounal  frequency,  S,  is  derived  in  a 
closed  analytical  form  as  follows: 


2V  o  r  /  r  \ 

s  -  V  °  M-V) 

"v»rc  \vr£/ 


(i.1) 


v/nere,  for  purposes  nerem,  we  estimate  the  characteristic  radial  dimension, 
r  ,  is  related  to  tne  waxe  width  (Refs.  6  and  43)  approximately  by, 


and 


U.2) 

Uo) 


(r.q) 


Por  vortex  snedding  m  tne  wake  of  a  cyliiider,  the  boundary  conditions,  along 
witn  the  Von  karman  spacing-ratio  (0.281),  enaole  determination  of  V  =  .77  V  » 
and  the  axi.al  growth  proportionality,  v  =  .36*  as  shown  on  Reference  1.  Tne 
variable-separable  form  of  solution  obtained  has  recently  been  verified,  during 
the  growth  period  (prior  to  vortex  breakdown),  at  least  partially  by  the 
analyses  in  Reference  13. 
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Use  of  a  Caucny  stretching  transformation  in  the  independent  variable,  r(t), 
jointly  witn  a  transformed  stream-function  P(0,r)  of  mixed  character  in  the 
orig_nal  dependent  and  independent  variables,  results  in  the  equation 


P  +  aiP2  +  o2P3  + 


+  f(t) 


(bo) 


Eor  radial  inflow  (see  Appendix  a),  which  is  characteristic  of  a  growing 
waxe  vortex,  R,.  is  negative.  Thus,  the  linearized  form  of  Equation  (Bo), 
dropping  ail  but  the  first  term  in  brackets,  results  in  eigenvalues  corresponding 
to  spatially  periodic  form  for  particular  roots  of  in  Equation  (E.3).  v/itu 
tne  hypothesis  of  characteristic  vortex  modes  (dipole,  quadropole,  etc. )f  the 
Reference  1  analysis  reduces  Equation  (B. 1)  to  a  form  for  S  pure_y  as  a 
lonotion  of  Reynolds  dumber  wmcn  agrees  quite  well  with  test  data  for  vortex 
shedding  behind  cylindrical  bodies. 


.oecause  of  the  nature  of  the  transformations  made  m  Reference  1 ,  the  solution 
.s  primarily  valid  in  outer  vortical  layers,  or  for  radii  beyond  tne  hard  core. 
It  is  ji  this  outer  region  wnere  tne  discussion  in  Section  3*2  of  this  Volume 
hypothesizes  that  R^  may  effectively  correspond  to  outflow,  1L(>0.  in  tuat 
case,  no  oscillatory  eigenvalues  exist  for  Equation  (Bo),  ana  characteristic 
solutions  w_il  tend  to  grow  without  bounds.  Conceptually,  then,  tms  could 
amount  to  a  necessary,  but  peiiiaps  not  sufficient,  condition  for  vortex  burst 


to  occur. 


